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1 Introduction 



1.1 Intention 

The aim of ^(oops is to calculate one-particle irreducible Feynman diagrams with one or two 
closed loops for arbitrary processes in the Standard model of particles and related theories. 
Up to now this aim is realized for all one-loop diagrams with at most three external lines 
and for two-loop diagrams with two external lines. 

The results can be returned both algebraically and numerically. All necessary tensor in- 
tegrals are treated for arbitrary masses and momenta. Those two-loop two-point functions 
for which no analytic result is known are integrated numerically by ^(oops. 

The calculations are performed in Maple V - a language for symbolic computations 
I2,|2|. Numerical integrations are done using VEGAS - a procedure for adaptive Monte 
Carlo integration ||§[]. For convenient input xloops has an Xwindows interface based on 
Tcl/Tk §]. 

The package consists of the following parts: 

t> Input via Xwindows interface: 

In a window the user selects the topology which shall be calculated. A Feynman 
diagram pops up in which the particle names have to be inserted. 

> Processing with Maple: 

The selected diagram is evaluated. The necessary steps for reducing the numerator 
(SU(N) algebra, Dirac matrices) are performed. The result is expressed in terms of 
one- and two-loop integrals. 

O Evaluation of one-loop integrals: 

One-loop one-, two- and three-point integrals are calculated analytically or numeri- 
cally to any tensor rank using Maple. This part was already released separately |],[7|. 

> Evaluation of two-loop integrals: 

All two-loop two-point topologies including tensor integrals are supported by the 
Maple routines. For those topologies where no analytic result is known, sloops creates 
either an analytic two-fold integral representation or integrates numerically with the 
help of VEGAS using C++ [p|. A parallelized implementation of the VEGAS algorithm 
invented by R. Kreckel ||[l(J is part of this distribution. 

The other topologies can be calculated analytically or numerically like one-loop 
integrals. 

In the second chapter we briefly denote how to install and start the program. Chapter 
three gives several examples for quick and easy usage. A more detailed description of the 
user-interface can be found in chapter four of this manual. 
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1 . 2 Copyright 

The copyright of this program package is owned by the Johannes Gutenberg-Universitat 
Mainz, Germany. 

7(loo-ps is " cite-ware" . That means that you can freely use and redistribute it, but you 
have to cite \Coops if you publish any results that you achieved with the help of it. 

Any further program development which intends to use the name sloops or the incor- 
poration of sloops or a part of the sloops code into any other program needs the permission 
of the authors. 

With the receipt of tfoop s the customer accepts that he uses the program entirely on 
his own risk. Neither the authors nor the distributors of \loops are liable for any loss or 
damage as a result of any use which is made of this program. 

1.3 Comments and bug reports 

sloops is still under development. So if you have questions, comments or if you have found 
any bug - please don't hesitate to contact 

xloopsOthep . physik . uni-mainz . de 

or one of the authors. 



2 Preliminaries 



2.1 Requirements 

^Coops can be used on any computer where Maple V is working. If Maple is not yet installed 
on your computer system you have to do that first. For numerical integration of two- loop 
integrals a compiler for C++ is necessary. The Xwindows interface needs Tcl/Tk and is 
at present only available for Unix. On other platforms or if Tcl/Tk is not installed the 
Xwindows interface will not work. In this case you can use the Maple part directly (cf. 



sect. and |2.5|) . The Maple part covers most features of the \Coops package but the input 



may be found not as user-friendly as if the interface is used. 



2.2 Installation 



The distribution of the tfoops package is accessible via WWW at 
http : //wwwthep . physik . uni-mainz . de/~xloops/ 
It consists of one packed (tared, gzipped) file 



xloops-1 . . tgz 

This file contains all files needed for sloops. To install the package unpack xloops-1 . . tgz: 



tar -xzvf xloops-1 . . tgz 



7(Coops then automatically creates a subdirectory of the current directory which is called 
xloops. The package is then deposited there. In addition there will occur several subdi- 
rectories of xloops. A complete list of all files in the distribution is given in appendix B.4 . 
Moreover, there exists an additional (tared, gzipped) file 



xloops-lib . tgz 



in the same directory as xloops-1 . . tgz which contains a library of one-loop diagrams 
(cf. sect. |4.5| ). For convenience one should unpack this file in the same directory where 
xloops-1 . .tgz was unpacked before, ^loops will read from this library whenever a one- 
or two-loop integral has to be evaluated. If this library is missing sloops will create the 
necessary files automatically, but this will usually take some time. If you have the bash 
shell type 



cd xloops/ 
configure 
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This will settle the configuration of the Xwindows interface, configure will request 

t> the release of Maple V (1, 2, 3 or 4) you are using. There were subtle changes from 
release 2 to 3 in handling local and global variables, therefore it is important to 
select the correct option. Since we tested only Maple V release 1 and 3, for release 2 
and 4 configure internally does the following replacement: 

release 2: use the same files as for release 1 

release 4: use the same files as for release 3 (with some patches) 

which is expected to be consistent. 

> the WWW browser which should be used for the on-line html version of this manual. 
The name of the executable (cf. netscape, mosaic, . . . ) must be inserted. 

D> the character encoding of your terminal (Isolatin 1 or standard ascii). 

[> whether multiprocessing is enabled. This question is necessary to choose between 
normal VEGAS or the parallelized version of VEGAS running on multiprocessor com- 
puters only. If multiprocessing is enabled the number of processors is requested. 

> the platform. Several frequently used systems (Linux, Digital Unix, SUN Solaris, IBM 
AIX) are distinguished from other Unix to install the right Makefile for the numerical 
integration. 

> whether a library of one-loop integrals should be build. It can take a lot of time 
(up to an hour) on slower computers to build the whole library. Therefore it is 
recommended not to do that, but to load a pre-built library xloops-lib . tgz from 
our WWW server instead. 

After that several checks are performed. First an embarassing bug occuring in Maple 
releases up to release 3 is tested: If a file is read in a procedure, no output of the procedure 
will appear. Therefore we have written a little workaround for this bug, which can be 
installed by the user. 

configure continues with a test run of the Makefile which is needed for numerical 
integrations. This Makefile depends on the platform. The run will take some time. 

Finally you are asked to agree that automatically a mail is sent to the developer team 
of the tfoops package. This mail is understood as registration of your xloo-ps installation. 



A typical run of configure is given in appendix A.l . configure automatically creates 
the sloops executable file. 

If you have root access you can make it executable for all users. You then can set a 
symbolic link to /usr/bin, i.e. 

cd /usr/bin 

In -s ( directory) /xloops/xloops . 

where {directory) is the directory where ^loops was installed. 

If you don't have bash you can find the necessary instructions in the file README which 
is part of the xloops distribution. 



2.3 Getting started 



You can now type 
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. /xloops 

from the directory where \loops was installed before. This starts your Troops session. The 
main window of xl°°v s (°f- fig- 3.1) will appear. If stoops was made executable for all users 



(cf. previous section) it can be started with 
xloops 

from any directory. 

There is on principle no restriction in the number of tfoops sessions which can be used 
simultaneously - by the same user or by different users. On the other hand the numerical 



integration with VEGAS should be used only once by each user (cf. sect. 4.6). 



2.4 Installing the Maple part only 

If you are only interested in the Maple routines you have to start the installation as in 



sect. 2.2, but you then skip configure and just copy all files with extension .ma 



> loops. ma 

> fmrules.ma 



> fmuser.ma 



> evalproc.ma 

> oneloop.ma 

> cfcn.ma 



> pv.ma 

> r.ma 

> simple. ma 

> twoloop.ma 

> numint.ma 

> values. ma 

> unvalue.ma 

> mess. ma 

in one directory instead and delete the rest of the distribution. The correct .ma files you 
can find in the subdirectory xloops/MapleVRl - if you use release 1 or 2 - or in the 
subdirectory xloops/MapleVR3 - if you have release 3 or higher. 
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2.5 Usage within a Maple session 

You can use \loops within an ordinary Maple session as well: After having started Maple 
just type the following lines: 

LoopPath : =(path) ; 
read' (path) loops .ma' ; 

(path) describes the directory where the .ma files are located. If the ^loops package is 
installed properly this should be the path of the ^(oops directory to which you have to add 
MapleVRl - if you are using release 1 or 2 of Maple V - or MapleVR3 - if you are using 
release 3 or higher. If the .ma files were copied to the current directory before you can skip 
the assignment of LoopPath and only type 

read 'loops. ma' ; 

Now you can use the functions described in chapter ||. 



3 Quick reference 



3.1 Main window 



After having typed the command xloops the main window (cf. fig. 3.1) appears on the 
screen. It consists of three parts: 

1. The menu bar which contains several menues for ^Coops commands (for details cf. 

sect. prop . 

2. The topology bar which displays all possible topologies for a definite m-loop n-point 
function. The selection of m closed loops and n external lines is done with the help 
of the menues Loops and Ext. Lines. 

3. The text window with scroll bar which contains the results of calculations. 



/ menu bar 



- topology bar 



File Model Loops /Ext. Lines Options 



Help 
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Maple V loaded ! 

> r e ad ( 1 /home /xlo ops/xloops /Map le VR3 /lo op s . ma 1 ) 
Warning: new definition for &* 



»»»»»»»»»»»»»»»»»» This is XLOOPS Version 1.0 

by Lars Brucher, Johannes Franzkowski, Dirk Kreimer 

©1993-1997 Johannes Gutenberg Universitat Mains, Germany 

bytes used=1040496, alloc=S51S12, time=0.43 
bytes used=2049840, alloc=13104S0, time =0. 94 



N Maple text window 



Figure 3.1: Main window 
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Before we will explain the functions and their meanings in detail in the next chapter, 
we demonstrate here as an example how to use xfoops for the calculation of a two-loop 
self-energy diagram and a one-loop vertex correction. 



3.2 Two-loop self-energy 

We now perform the calculation of a diagram which contributes to the Higgs self-energy £ 
at the two-loop level: 




The underlying topology is the master topology, which is located in the fourth place in 
the topology bar (cf. fig. [O]). In this case the number of external lines and the number of 
closed loops is already fixed correctly, so one now has just to click on the topology. 

Then the master topology appears in a separate window, the diagram window (cf. 
fig. [T^). In this window for each internal and external particle the corresponding particle 
name has to be inserted. Only particles can occur here, no antiparticles. The allowed terms 
for particle names are listed in sect. 4.2.4[ They can be obtained alternatively with the 
menu Help. The flow of charges and fermion or ghost numbers is determined by the arrows 
of the lines. The direction of an arrow is reverted with a click on it. 



Sub topologies 



no 



Export 




Evaluate I Eval. Full Eval. Numeric Close 



Figure 3.2: The diagram for the two-loop Higgs self-energy 
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Since in our example we have to deal only with chargeless scalars, the direction of 
arrows is meaningless and we have to insert the name H or higgs for each particle only. 

For the sake of clearness, the flow of external momenta in our conventions is displayed 
in blue. In this example there is only one momentum 

91 = (910,0,0,0) = (<? 10 ,0). 

Vanishing components are put together to a common zero-vector. In our case there exists 
only one parallel space component qiQ. The orthogonal space is represented through a 
common 0. 

When all particles are inserted, Troops is ready for calculation. The diagram window 
gives you the possibility to select different modes of evaluation. 

> Click on the ( Evaluate | button: The output - for simplicity - is returned in terms of 
Oneloop and Twoloop functions which correspond directly to one-loop and two-loop 



functions explained in sect. \LQ and |4.4| . In our example the result consists only of 



the scalar master two-loop two-point function which %Coops calls Twoloop2Pt2: 

> Gl := EvalGraph2 (7, [higgs , higgs , higgs , higgs , higgs , higgs , higgs , higgs , higgs , 

> higgs , higgs , higgs] ) ; 

Gl := [ 

81 84 2 2 2 222 

CI = I Mhiggs e (1+2 eps Ln(4 Pi MU ) + 2 eps Ln(4 Pi MU ) ) 

8192 

Twoloop2Pt2(0, 0, 0, 0, 0, qlO, Mhiggs, Mhiggs, Mhiggs, Mhiggs, Mhiggs) 

/ 4 4 8 

/ (sin(tw) Mw Pi ) , 

/ 

CI] 

Generally, the output is decomposed in one or several form factors CI, C2, ... to 
reflect the Dirac 7 or Lorentz structure of the result. Therefore the result is a list of 
form factors. The last entry in this list is the defining equation for the form factors. 
In our example of a scalar self-energy there appears only one form factor CI and the 
defining equation is just the scalar CI. A more complicated situation can be found 
in the next section (for details cf. sect 4.2.10| ). 



Moreover sloops assigned the whole result to a variable Gl. Within an sloops session 
each calculation will get such a name in an unique way. This gives the user the 
opportunity to perform further calculations with the results of diagrams, for instance 
to sum up all evaluated graphs which contribute to the same process. 



> Click on the [ Eval. Full] button: ^(oops evaluates directly the Oneloop and Twoloop 
integrals. Now each form factor CI, ... is a Laurent expansion in terms of the ultra- 
violet regulator e, where the significant coefficients of this expansion - C(e _1 ), O(e ) 
in the one-loop case, 0(e~ 2 ), 0(e~ 1 ), O(e ) in the two-loop case - are denoted in 
form of a list (cf. sect |4.2.10| ) . 

For those two-loop topologies where no analytic result is known - like in the case 
of our example, the finite part contributes two entries to the output list. The first 
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of them is a list itself which denotes an integral representation, the second contains 
the part which is analytically calculable: 

GFl := [ 
CI = [0, 0, 

8 4 

81 Mhiggs e 
[ > 

8192 4 4 4 

sin(tw) Mw Pi 

2 2 

2 Mhiggs 2 Mhiggs 

I Ln(Sqrt(x + I rho) + 7.1 + Sqrt(y - 2 y + 1 + I rho)) 

2 2 

qlO qlO 
4 

2 

qlO (- 2 x - 1) (- 2 y + 1) 
2 

2 Mhiggs 

- 4 I Ln(Sqrt(x + 2 x + 1 + I rho) + '/.I 

2 

qlO 

2 

2 Mhiggs / 2 

+ Sqrt(y - 2 y + 1 + I rho)) / (qlO (- 2 x - 1) (- 2 y + 1)) 

2 / 

qlO 

2 2 
2 Mhiggs 2 Mhiggs 

I Ln(Sqrt(x + I rho) + + Sqrt(y + I rho)) 

2 2 
qlO qlO 
_ 4 + 4 

2 

qlO (- 2 x - 1) (- 2 y + 1) 

2 2 
2 Mhiggs 2 Mhiggs 

I Ln(Sqrt(x + 2 x + 1 + I rho) + + Sqrt(y + I rho)) 

2 2 
qlO qlO 

2 

qlO (- 2 x - 1) (- 2 y + 1) 



x = - infinity . . infinity, y = - infinity . . infinity] , 

0], 

CI] 
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2 

2 2 Mhiggs 

11 := Sqrt(x + 2 x y + y + I rho) 

2 

qlO 



In our example the result is finite, therefore the first two entries in the list are zero. 
The only nonvanishing contribution comes from the integral representation. A more 



detailed description of the integral representation is given in sect. 4.4.3 . 

The result is a function of the Higgs mass Mhiggs and the external momentum qlO. 
The coupling of the vertices involve the W mass Mw, the Weinberg angle tw and the 
electromagnetic coupling constant e. rho denotes the infinitesimal imaginary part 
of the masses in the denominator (cf. sect. 4.2.1C| ), x and y are the variables which 
have to be integrated out numerically. 



> Click on the [ Evai. Numeric] button: The integral representation described above is 
evaluated directly with the help of VEGAS. This option exists only for two-loop 
integrals, because all one-loop integrals are analytically calculable. 

Please be aware that ?cCoops can only perform the numerical integration if all variables 
- Mhiggs and qlO in our example - are assigned with numerical values. For that 
purpose the menu File provides the Insert Maple Command entry. If this option is 



selected, a window (cf. fig. ^3) appears, where any command can be passed directly 
to Maple. In our example we insert 



ql0:=91. 17; 
Mhiggs : =300 ; 

to evaluate the integral on the Z resonance for an assumed value of 300 GeV for the 
Higgs mass. After clicking on I Evai. Numeric I ^oops starts the numerical integration. 







Evaluate Maple Command 


FjI 


Please insert your Maple command: 


|q10 := 91.17; Mhiggs := 300; 


Evaluate 
U 




□ear Buffer I 


Close I 

- 



Figure 3.3: Direct Maple input of masses and momenta 

The result including the numerical uncertainty appears in the main window and 
reads: 

GMl := [ 
ci = [o, o, 



16 e 



[.66598680119771701375*10 



4 4 
sin(tw) Mw 
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-8 -5 -7 -7 

[ - .354349*10 + .880157*10 I, .13042*10 + .148817*10 I]], 0] , 



CI] 



In the result - compared to the ( Evai. Full | option - the integral representation is 
now replaced by the numerical result of the integral. This result consists of a list: 
a pre-factor and another list. This list contains two entries: first the result of the 
VEGAS calculation, then the uncertainty VEGAS detected for this result (for details 
cf. sect. plO| ). 

In our example the real part of the diagram is less than the numerical uncertainty, 
whereas the imaginary part is significantly larger. This means that the real part in 
fact is zero. Since we calculated the entire diagram which is not precisely E but — iE, 
this is what one expects below threshold. 

Alternatively to the direct input of Maple commands described above it is possible 
to assign all known masses of particles, couplings and elements of the CKM matrix 
with the option Insert Particle Properties of the Options menu. 

The input information for VEGAS can also be determined with the menu Options: 
Click on Numeric and an additional window will appear where the input parameters 
for VEGAS can be fixed (for details cf. sect. 4,1.4| ). 



Finally the result can be saved in text mode to a file. Use in the menu File the entry 
Save Maple Output. 



\> Click on the [ Close] button: The diagram window disappears. 

The four buttons described above are completely independent. You can for example select 
( Evai. Numeric ) without having clicked on [ Evaluate | or [ Evai. Full] before. 



3.3 One-loop vertex correction 

What follows is the calculation of a one-loop contribution of the correction to the eeZ 
vertex: 




First one has to select the correct values in the menues Loops and Ext. Lines of the main 
window to get access to the one-loop three-point topologies. A click on the correct topology 
opens the diagram window (cf. fig. 3A) where the particle names have to be inserted. Again 
the correct particle names can be found in sect. 4.2.4 or in the menu Help. The direction of 
the fermion line must be adjusted by clicking on the arrows. The flow of external momenta 
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I 



Subtapologies 



Export 



^ 

iort 



(g2Q,g2i^ 




(q10,0). 



(q10-qg0,-q21 ] 



Evaluate 



E. Full E. More Close 



Figure 3.4: The diagram for the eeZ one-loop vertex correction 



is displayed in blue again. The three momenta are 

91 = (910,0,0,0) = ( 9 io,0) 

92 = (920,921,0,0) = (920,^21) 

9i - 92 = (910 - 920,-921,0,0) = (910 - 920,-921) 

where the vanishing orthogonal space components are not displayed in the window. 
910,920 , 921 are the kinematical variables tfoops uses for three-point functions. To get a 
numerical result, we assign these variables with 

ql0:=Mz0; q20:=MzO/2; q21 : =evalf (sqrt (Mz0"2/4-Melec"2) ) ; 

using the option Insert Maple Command of the menu File so that all external particles 
are on-shell. Moreover the masses and couplings are assigned by the option Insert Particle 
Properties of the Options menu. The renormalization scale MU and the unit GeV can be 
omitted by setting 

MU:=1; GeV:=l; 



Clicking on the [ E. Full j button gives the following result: 

GF2 := [ 

CI = [.000010058632263504308717 I, 

- .00086276496443282788638 - .0072236143340212361270 I], 
C2 = [.00010510130634321490041 I, 

- .0090149159887518402249 - .075478582279918574719 I], 
C3 = [0, .00083116483915048522171 + .0071524795241752712917 I] 
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C4 = [0, .0086847304974275009526 + .074735304153769864092 I], 

-9 -8 
C5 = [0, .40863107990668029069*10 + .35164209243516742716*10 I], 

C6 = [0, .000018229897665296319548 + .00015687498271956514041 I], 

C7 = [0, .00019048176819004810260 + .0016391657617525244669 I], 

-28 -27 
C8 = [0, - .25407817700937862122*10 - .21701516787925566779*10 I], 

C9 = [0, .000017536874322335126422 + .00015153965406412662322 I], 

CIO = [0, .00018324045975332954356 + .0015834176245538375947 I], 

-9 -8 
Cll = [0, - .20431553995334014636*10 - .17582104621758371369*10 I], 

-8 -7 

C12 = [0, .22971899930888218199*10 + .19620932339352982103*10 I], 

-25 

C13 = [ - .24074590585723032270*10 I, 

-6 -5 
.34651167097341202104*10 + .28882798115295165844*10 I], 

-5 

C14 = [0, .36206542130597753444*10 + .000030179250351750043557 I], 

-28 -27 
C15 = [0, .12703908850468931061*10 + .10850758393962783388*10 I], 

-6 -5 
C16 = [0, .34651167090811890878*10 + .26676643227928772485*10 I], 

-5 

C17 = [0, .36206542123775363980*10 + .000027874068547868333215 I], 

CI (1 &* Dg(nul)) + C2 &*(1, Dg5, Dg(nul)) + C3 &*(1, DgO, Dgl , Dg(iml)) 

+ C4 &*(1, Dg5, DgO, Dgl, Dg(nul)) + C5 (1 ft* ONE) q2(nul) 

+ C6 (1 ft* DgO) q2(nul) + C7 &*(1, Dg5, DgO) q2(nul) 

+ C8 &*(1, DgO, Dgl) q2(nul) + C9 (1 ft* Dgl) q2(nul) 

+ CIO &*(1, Dg5, Dgl) q2(nul) + Cll (1 ft* ONE) ql(nul) 

+ C12 (1 &* Dg5) ql(nul) + C13 (1 &* DgO) ql(iml) 

+ C14 &*(1, Dg5, DgO) ql(nul) + C15 &*(1, DgO, Dgl) ql(mil) 

+ C16 (1 ft* Dgl) ql(nul) + C17 &*(1, Dg5, Dgl) ql(nul) ] 

The output is decomposed in several form factors CI, C2, ... to reflect the Dirac 7 and 
Lorentz structure of the result. The last entry in this list is the defining equation for the 



form factors (cf. sect. 4.2.10) 



4 Detailed description 



There are two intentions covered by xioops: 

> Either the calculation of complete Feynman diagrams. Here you need the Evaluate 
Graph window (see sect. |4.1|) where you can insert particles into a Feynman diagram. 



Alternatively you can use the EvalGraph procedures in Maple described in sect. [4J2 



D> Or the solution of single integrals occurring in the calculation of one- and two-loop 
diagrams. Then you can make use of the OneLoop and TwoLoop procedures explained 



in sect. n.3i and 4.4. 



The creation of a library for the integrals follows in sect. kL3. The numerical integration 



of two-loop integrals is explained in sect. 4.6 



4.1 The Xwindows interface 
4.1.1 The functions of the menu bar 

The menu bar provides a variety of different functions for the user. These functions are 
combined to groups which appear as menues if the user clicks on it. The list below gives 
a short overview. Additional descriptions can be found in the following sections. 

> The menu File (cf. fig. f4.1| a) has the following entries: 

> Insert Maple Command allows the insertion of an arbitrary command which is 



passed directly to Maple (cf. fig. |3.3|) . Especially the numerical input of masses, 
momenta and renormalization conditions can be done easily this way. 

> Save Maple Output saves the output as it is denoted in the Maple text window to 
an extra file. The name of the file has to be determined in a separate window, 
which appears if the user clicks on this entry. 

> New Process allows the user to collect a set of graphs and save them with a 



common process name. Detailed explanation can be found in section 4.1.3. 

> With Load Process processes which were declared with New Process before can 
be reloaded into the current Troops session (cf. sect. |4.1.3| ). 

> Create Process will call the graph generator to produce all graphs to a given 
number of loops and external lines in a pre-defined model (in future versions of 
Troops only). 

> Restart Maple terminates the current Maple session and removes all assigned 
variables. Maple and tfoops are started again. 

> Quit terminates xfoops. 
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Insert Maple Command 
Save Maple Output 



New Process 
Load Process 



Create Process 



Restart Maple 
Quit 



(a) 



About 

Help 

Tutorial 



(b) 



V History 

W Oneloop Library 



V Insert Particle Properties 
■ Set Neutrino Masses 
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Figure 4.1: The menues of the menu bar 
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\> The menu Model (cf. fig. [4.l| d) offers several physical models for calculation (for 
details cf. sect. 4.2.3| ) . 



l> The menu Loops (cf. fig. |4.l| e) selects the number of closed loops for the calculation. 
The user can choose between 0, 1 and 2. A change of this value results in different 
topologies displayed in the topology bar. 

t> The menu Ext. Lines (cf. fig. [4.l| f) allows the user to choose the number of external 
lines. At the moment the upper limit is three. The change of this number results in 
a different topology bar as well. 

[> The menu Options (cf. fig. [4.1| c) contains the following entries: 

t> If History is selected tcCoojjs saves all calculations under a common process name. 
Those diagrams which were already saved before are read from disk. Otherwise, 
if History is not selected they are recalculated (cf. sect. 4.1.3| ). 



t> With the selection of Oneloop Library ^[oops takes all one-loop functions from 
a pre-built library. Otherwise the one-loop functions are calculated every time 



they are needed (cf. sect. 4.5 and appendix 

If Insert Particle Properties is selected the values of masses and couplings reported 
by the Particle Data Group |11] are inserted (cf. sect. [4.2.4 ). 



> Set Neutrino Masses assigns all neutrino masses to zero (default) or, if not se- 
lected, leaves the neutrino masses unassigned. In this case an undefined neutrino 
mixing matrix exists as well (cf. sect. 4.2.4] ) . 
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> Massless Light Fermions assigns the masses of all fermions except the b and t 
quark to zero (cf. sect. 4.2.4| ). 



> Numeric defines the number of points and the number of iterations for the 



numerical integration with VEGAS (cf. sect. 4.6). 

> Window Options allows the user to adjust the appearance of the Xwindows in- 
terface, especially the placement of the red arrows of the diagram window (cf. 
sect. |4.1.4j ). 



> Save Options saves all options in the file .xloops in the home directory of the 
user. 

t> The menu Help (cf. fig. |4.1| b) has three entries: 

> About shows a brief information on tfoops. 

> Help gives a short help where especially the notation of all particle names in 
xioops is written. 

> Tutorial calls a WWW browser which displays the on-line html version of this 
manual. 

4.1.2 Topology bar and diagram window 

According to the pre-defined values for the numbers of closed loops and external lines the 
topology bar displays all possible topologies. If the user clicks on one of them this topology 
appears in a separate window, the diagram window. This window consists of three parts: a 
menu bar, an area for inserting particles (diagram area) and a command bar with different 
buttons for evaluation (cf. fig. |4.2|) . 

The menu bar provides two options. The right one, Export, gives the user the possibility 
to export the graph to a Postscript file. Only the diagram area of the diagram window is 
exported, menu bar and command bar are not included. 

The left option, Subtopologies, will only be available in future versions. It will provide 
the possibility to switch between several subtopologies, which are obtained when external 



lines are interchanged. The upper topology in fig. 4.2 may serve as an example. This 
topology can be rotated by 120° or 240° which corresponds to a cyclic permutation of the 
external lines. This changes the topology to the same type of topology but with different 



momentum flow. In contrast, the lower topology in fig. |4.2| does not change if it is rotated 
by 120° or 240°. So this topology does not exhibit different subtopologies. 

In the diagram area the topology is displayed graphically. For every propagator an 
entry for the insertion of particles is located. The correct terms for particle names are 
written in section 4.2.4J . They can also be obtained from the Help function of the 



mam 



window. The flow of charge and fermion and ghost numbers are determined by the red 
arrows on the propagator lines. Clicking on an arrow changes its direction. The correct 
flow of charge and fermion and ghost numbers can be settled in this way. Only particles, 
no antiparticles can be inserted. Charge conjugation may be achieved by inversion of 
arrows. Please keep in mind that the arrows don't describe the direction of momenta. The 
convention of momentum flow is shown graphically in the diagram window in blue. It is 
defined for all diagrams in the same way, to facilitate a final summation of all calculated 
diagrams. 

If the user has inserted all particles and checked the charge flow, the command area 
gives him several possibilities for the evaluation of the diagram. They correspond to the 
different modes of evaluation as explained in section 4.2.10| . 
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Figure 4.2: Different diagram windows 
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4.1.3 Loading and saving of processes 

For the evaluation of complete processes it is possible to save the contributing graphs with 
a common process name. Therefore Tcfoops provides the File menu entries Load Process and 
New Process. 

If Load Process is selected the file window (cf. fig. |4.3j ) appears on the screen. In this 





"-I 


Existing Saved Processes 


Higgsjilus 
qed 


H+ -> W+ hO 

The basic OED graphs 


A 


Process Name: 


Iqed 








OK | 


Cancel | 









Figure 4.3: The Load Process file window 



window the processes which were already saved are listed on the left hand side. In addition 
for every process a user specified description is given on the right hand side. This may 
help to reidentify a process later. To select a specific process one clicks on a process name 
on the left. After clicking on the [ OK | button the history table (cf. fig |4.4|) of the selected 
process appears whereas ( Cancel | closes the file window. In the history table all diagrams 
which were included in the corresponding process occur with a button. A click on such 
a button opens a diagram window (cf. fig p~^ ) with which the Feynman diagram can be 
evaluated. The history table contains two menues: 

t> In Graphs one has the opportunity - if the number of diagrams exceeds the capacity 
of the window - to switch to another sheet of diagrams. The same effect have the 
arrows at the bottom of the history table. 

> Close lets the history table disappear. 

If the entry History of the Options menu is selected all results are saved. With History 
all calculations which were already saved in this or another ^Coops session are read in 
automatically and not calculated again. For each process xioops automatically creates a 
subdirectory of the directory xloops_user where the results are saved to different files. 
If History is selected without having specified a process before ^loo-ps assumes that a 



new process shall be added to the list and opens a similar file window (cf. fig 4J3) - the 
same window appears if New Process is selected. Two entries are provided for inserting a 
new process name and description at the bottom of the window. After the user inserted the 
name and description for his process and clicked on the [ OK ] button a new - and therefore 
empty - history table appears on the screen. All diagrams which are calculated from now 
are included in the history table as additional entries. The topology of the graph and the 
list of interacting particles are displayed. 
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Figure 4.4: The history table containing diagrams of the process H — > H 
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Figure 4.5: The New Process file window 



4.1.4 Window options 



With the entry Window Options of the Options menu the appearance of some windows of 
xCoops can be adjusted to fit the terminal characteristics (cf. fig. |4.6| ). First, the size of the 
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Set arrow place 



Set Factor 



Close 



Figure 4.6: Window options of tfoops 



history table (cf. fig. 4.4) can be changed by decreasing or increasing the numbers of the 
Rows and Columns entry in fig. 4.6. Then, the position of the arrows of the propagators in 
the diagram window can be corrected. This might be necessary on some terminal types, 
where the arrows of the propagators are misplaced. This has to do with the fact that the 
relation between screen dots (as used by Tk) and centimeters is different on different X 
terminals. The embarrassing consequence can be that the arrows of the propagators are 
misplaced. The problem can be avoided by inserting a global Scale Factor. Alternatively, 
if this factor is not obvious, \loops will calculate it. Therefore the user must insert the 
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position of the arrow on the ruler in the Arrow Scaling field of fig. 4.6.^] 

4.2 Evaluation of Feynman diagrams 
4.2.1 Input 

To evaluate a particular Feynman diagram itfoops needs to know four ingredients: 

1. The number of closed loops. 

2. The number of external lines. 

3. The topology, that means the information how the lines of the diagram are connected. 

4. The particles of the diagram. 

The first two numbers have to be selected with the help of the menu bar entries Loops and 
Ext. Lines (cf. sect. |4.1.1|) . Then the topology can be selected on the topology bar so that 
the diagram window pops up where the particles can be inserted (cf. sect. 4.1. 2| ) . 



Internally the input from the Xwindows interface is converted to a procedure which 
performs the calculations in Maple. It is also possible to start directly from this routine 
for Maple - either with the Insert Maple Command of the File menu or within an ordinary 
Maple session where the Maple part of xloops was read in before (cf. sect. |2.5| ). 

For that purpose the internal Maple commands will be mentioned in the following. All 
commands which follow without reference to the Xwindows interface have to be passed to 
Maple. 

The input is then of the following form 

> for tree diagrams: EvalGraphO (n , (list)) ; 

> for one-loop diagrams: EvalGraphl (n, (list)) ; 

> for two-loop diagrams: EvalGraph2(n, (list)) ; 

The different topologies with the same number of loops are numbered. The number of the 
desired topology has to be inserted in the function call for n. The correct number for each 



topology will be given below in sect. 4.2.2| . (list) describes the list of particles involved in 



the Feynman diagram. The conventions for the ordering of the particles are also explained 



for every topology in 4.2.2 



4.2.2 Internal Notation 

What follows is a list of all topologies which at present can be solved by ^Coops. This list 
is only of interest for direct Maple input. To fix the flow of charges and of fermion or 
ghost numbers in an unique way, each internal propagator is assigned with two particle 
names. For example [up.upbar] and [upbar.up] describe the up-quark propagator but 
with different direction of the quark line arrow. External lines are described by only one 
particle name. 

As a convention all particles at each vertex are incoming. The arrows of the diagrams 
displayed below only describe the flow of momenta. Of course these arrows need not to 
coincide with the arrows of fermions or other charged particles. 

For each topology we give the correct EvalGraph command as well as the command 
for the corresponding OneLoop or TwoLoop integral. 

1 If the scaling factor is already correct, the arrow should point on the 5. 
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Tree level 

Three-point diagrams 

EvalGraphO(l, [®,©,®]) 




Four-point diagrams 



EvalGraphO(2, [®,©,®,@]) 
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One-loop level 
One-point diagrams 



EvalGraphl(l, [©,©,©]) 




OneLooplPt(p, to) 



Two-point diagrams 

I 1 



EvalGraphl(2, [©,©,©,©]) 
I 




OneLooplPt(p, to) 



4 Detailed description 



29 



EvalGraphl(3, [©,©,©,©,©,©]) 




I 






/© 


©\ 


9i I 

> ♦ 


\ 9i 


© \ 


/ © 

©/ 
mi ^/ 

Z + 9i 


0neLoop2Pt(p ,P±, Qio, m 1 , m 2 ) 



Three-point diagrams 



EvalGraphl(4, [©,©,©,©,© 


,©,©]) 


I 








/© ©\ 


© /< 12 


9i / \ 




> f 




© \ y 




V© ©/ 


(2)\Q2 - qi 


"21 








DneLoop2Pt(p ,P±, 9io, mi 


m 2 ) 



EvalGraphl(5, [®,(| 


),©,©,©, 


© 


,©,©,©]) 








-®» 




©^ 








m 3 © 










\l 


+ 92 


^ © 






i + 


^^^^ 




92 - 9i 








< 








© 


0neLoop3Pt (po, Pi , P± 


, 910, 920, 921, 


mi,m 2 ,m 3 ) 



30 



4 Detailed description 



Two-loop level 
Two-point diagrams 





EvalGraph2(6, [®,( 


2), 




),©,®,®,®,®,©D 










k 


















I 

Qi I 
> i 


/©^ 

® 


© 


® 

I + k 


^®\ i 

®\ 91 




© 1 






mi ^ 
l + Ql 


f > 

® / 




TwoLoop2Ptl(p ,P± 




,r±,s, qio, mi,ni2,ni3, 777.4, ^5) 





EvalGraph2(7,[®,©,®,@,©,©,©,®,®,®,®,©]) 




TwoLoop2Pt2(p , p±, r , r±, s, q w , mi, m 2 , m 3 , m 4 , m 5 ) 



4 Detailed description 



31 



EvalGraph2(8, [©,©,©,© 


©,©,©,®,®,®D 


a»3 y® 




/ "1,4 ✓ 

/© >fc 

yi i 


©\ 

1 "1 




©/ 






J + 91 




TwoLoop2Pt3(pO)P±> ^0) r ±) 


gio,mi,m 2 ,m3,m 4 ) 



EvalGraph2(9, [©,©,©,©,© 


©,©,©]) 




i + k 






m 2 v. 




/© ^ 
91 / © 


^m 3 © 


®\ 

\ 9i 




k 

mi 
l + qi 


® / 


TwoLoop2Pt4(p ,P± 


,r ,r±,s,q w , 


mi,m 2 ,m 3 ) 



32 



4 Detailed description 



EvalGraph2(lO,[®,©,®,@,© J ©,©,®,®,®]) 
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For each encircled number one particle has to be inserted. The correct symbols for all 
allowed particles are listed in the next section. The conventions for the OneLoop and 



TwoLoop functions are denoted in sect. 4.3 and 4.4. 



4.2.3 Models 

?c(oops distinguishes several models. The menu Model of the menu bar serves for the selection 
of the desired model. Internally the global variable Model selects the model which shall be 
used. The following models are declared in \loo-ps: 
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> Standard Model - the (minimal) standard model (electroweak and QCD). This cor- 
responds to Model :=SM; and is the default value 

> 2 Higgs Doublets Model - the standard model with two Higgs doublets which corre- 
sponds to five physical Higgs particles; Model :=THDM; 

t> QCD - QCD with three generations of Quarks; Model :=QCD; 

O Electroweak Theory - the electroweak sector of the standard model; Model :=EW; 

> QED - QED with three generations of Leptons and Quarks; Model :=QED; 
D> Phi"4 Theory - <^> 4 -theory in four dimensions; Model :=phif our; 

> Phi"3 Theory - (/> 3 -theory in four dimensions; Model : =phithree ; 

D> Goldstone Boson Equiv. Model - the equivalence model for Goldstone bosons - 
standard model vector bosons are replaced by scalar Goldstone modes (l^ O ; 
Model :=GBEM; 

> a user-defined model where arbitrary Feynman rules can be declared. To work with 
this option the user has to set Model :=user; and to read an additional file 

read' (path)imuser .ma' ; 

where {path} is the same path used for the assignment of LoopPath at the beginning 
of the Maple session (cf. sect. |2.5[ ), In the file fmuser.ma the user-defined rules have 
to be inserted. Initially all QCD Feynman rules are denoted there as an example. 
The conventions of declaration can be taken from this example. 



4.2.4 Particles and Masses 

What follows is a list which contains all particles Troops knows. The names valid in the 
Xwindows interface are denoted (second column from the left) as well as the names for 
internal use within Maple (third column from the left). In addition for each particle the 
models in which it is declared (right column) and the name which %Coops uses for the 
corresponding particle mass (second column from the right) are listed. 
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Please note that the masses of the neutrinos are set to zero by default. To keep them as 
variables, one can switch off the entry Set Neutrino Mass of the Options menu or use 

UnsetNeutrinoMassOO ; 

within Maple which removes any assignment of Mnuelec, Mnumu and Mnutau. As a conse- 
quence now a leptonic CKM matrix exists which behaves analogous to the usual CKM 



matrix of quarks (cf. sect. 4.2.6 ). The reverse command 
SetNeutrinoMassOO ; 

or the reselection of Set Neutrino Mass sets all neutrino masses back to zero. The leptonic 
CKM matrix equals now the unit matrix. 

The masses of photon and gluon are also set to zero by default. To keep them as vari- 
ables you have to type 



Mgamma : = ' Mgamma ' ; Mgluon : = ' Mgluon ' ; 

The masses of leptons and light quarks - all standard model fermions except b and t quark 
- can collectively be set to zero with the entry Massless Light Fermions of the Options menu 
or with 



SetLightMassOO ; 
and be unassigned with 
UnsetLightMassOQ ; 



or another click on Massless Light Fermions. Apart from this, masses, momenta and cou- 
plings - as they are explained in the following sections - can be assigned by the user. 

Moreover, if you want to assign all particle massesf] as well as the couplings, the 
Cabibbo-Kobayashi-Maskawa matrix (CKM matrix) and the Weinberg angle with their 
experimentally determined values, you can use the entry Insert Particle Properties of the 
Options menu or type 

read ' (path)va.lues .ma'; 

where {path) is the same path used for the assignment of LoopPath at the beginning of 



the Maple session (cf. sect. 2.5). The file values. ma contains the parameters according to 
the currently published data of the Particle Data Group [P-ljj. To remove all information 
of values. ma in a Maple session you can type 

read ' (path) unvalue .ma'; 

or switch off the Insert Particle Properties entry. Instead of the fermion names which were 
already introduced one can also use the more general terms 



upl, up2, up3, downl, down2, down3 . 



The difference is, that for these terms no electric charge Q and no weak isospin T31 is 



2 Unfortunately the Higgs mass is not yet included 
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defined. While \loops assumes 

Q(elec)=-1; T31(elec) = -± ; Q(up) = |; T31(up) = |; 

it has no values for Q(upl), T31(upl), . . . and the user has the opportunity to make use 
of this freedom. 

4.2.5 Feynman rules 

Since the conventions for Feynman rules vary in the literature, it is important to explain 
one's own rules: We use Feynman gauge, for electroweak processes according to the rules 
given in fllEfl , but with all three- vertices and all propagators multiplied by —i and all four- 
vertices and the integration measure multiplied by i. This convention coincides with the 
rules of - as far as these are applicable in Feynman gauge. From [16] we took the rules 



for QCD and the notation of the CKM matrix, too. The rules for the Two Higgs doublet 
model are from fl7| . 

Anyway, you receive the rule for one particular vertex easily if you just use the tree 
level topologies or call the EvalGraphO procedure. 

In our convention the electromagnetic coupling e and the strong coupling gs are posi- 
tive (the electric charge of the electron is -e). xfoops automatically inserts a factor (27r) 2e ~ 4 
for every closed loop and an additional — 1 for each fermion and ghost loop. This means 
that the dimension - which tfoops calls D - equals 4 - 2 eps. The correct symmetry factor 
of the diagram is also included. The scale which keeps the coupling constants dimensionless 
is called MU. 

4.2.6 Mixing angles and CKM matrix 

Three mixing angles are known by ^Coops: The usual Weinberg angle tw and - only relevant 



in the case of the Two Higgs doublet model - the angles alpha and beta (cf. [17Q. The 
elements of the CKM matrix are denoted by CKM (up, down) . . . , its complex conjugates 
are CKMC (up , down) .... If neutrino masses are not set equal to zero there exists also a 
leptonic CKM matrix with elements like CKM(nuelec , elec) . If the neutrino masses are 



set to zero - which is the default (cf. sect. 4.2.4 ) - this matrix equals the unit matrix. 



4.2.7 Momenta and Metric 

The external momenta ql(nul) , q2(nul) , . . . may occur in two different contexts, 

1. as Lorentz vectors ql(nul), . . . with Lorentz indices nul, . . . , 

2. contracted with another Lorentz vector. 

In the first case the momentum will appear in the output as it is. In the second case the 
contraction will be written in terms of parallel space components: 



nin ..< -. „ on 9i • 92 . / (gi • Q2? 2 
qlO=Vyi; q20 = — =-• q21 = W = q 2 ; 



(4.1) 



„™ qi ' qs ■ SllSi- / (gi • q^) 2 (<?2 • g 3 ) 2 2 

q30 = — q31 = -== ; q32 = W § 1 ^ <?3 

Vyi V-Q2 V 9i 92 



where the auxiliary property 



9i ' 92 f . , 

92 = Q2 — 9i (4.2) 

9i 
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was used. In appendix |X| we demonstrate how this works in practice. 

The metric tensor of nloops is G(nul,nu2). It obeys the usual conventions 

G(0,0)=1; G(l,l)=-1 

The contraction of G(nul ,nu2) depends on whether the indices nul and nu2 are orthogonal 
space indices - then G (nul , nul) = D — DimP - or belong to the entire Minkowski space - 
this means that G(nul,nul) = D. 

The totally antisymmetric tensor is Eps(nul,nu2,nu3,nu4) with Eps(0, 1,2,3) = 1. 

4.2.8 Dirac matrices 

Dirac 7-matrices ^loops calls Dg(mul), .... Like the external momenta they can appear 
either as Lorentz vector as they are, or contracted by another Lorentz vector. In the latter 
case they are denoted by 

DgO=^; Dgl = — iLj; Dg2 = - (4.3) 
V9i V-<?2 V-<?3 

with 



9l -93 92 • 93 . . \ 

93 = Qs — 9i 92 (4.4) 

9i 92 

and the abbreviation ( |4.2| ). The Dirac structure is written as a string of the non- 
commutative product &*. Since the program is also designed for diagrams with more 
than one fermion line it labels each string of Dirac matrices with a number. For example 
in the case of fermionic self-energies 

ft* (1, ONE); &*(l,DgO); &*(l,Dg5); &*(l,Dg5,DgO) 

will occur in the output - the latter two of course only in the case of W, Z or Goldstone 
boson exchange. ONE represents the identity of the Clifford algebra, Dg5 is 75. The first 
entry in the non-commutative string gives the number of the fermion line. The 1 means 
that the string belongs to the first fermion line which in this case is the only fermion line. 
In appendix |A] this treatment can be seen in practice. 

4.2.9 SU(N C ) algebra 

In QCD calculations ?cfoops keeps the colour number Nc unassigned by default, but with 
Nc:=3; 

the number is fixed to its standard model value. The SU(N C ) algebra needs as ingredients 

t> the identity matrix of the SU(N C ) algebra delta3(tal ,ta2) which should not be 
mixed with Kronecker delta8(tal ,ta2) which contracts the Gell-Mann matrices 
with the labels tal and ta2. Of course delta3(tal ,ta2) and delta8(tal ,ta2) are 
both diagonal with entries 1 in the diagonal. This means that 

^delta3(tal,tal) = N c = 3 ; V] delta8(tal ,tal) = iV 6 2 - 1 = 8 . 

tal tal 
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> the Gell-Mann matrices Ts(tal,ta2,ta3) themselves, where tal - which is the 
number of the matrix - can take values from 1 to iV| — 1. ta2 and ta3 are the 
indices of a particular Gell-Mann matrix and therefore have values from 1 to N c . 
The Ts(tal ,ta2,ta3) are normalized by the following trace condition: 

Ts(tal,ta3,ta4) Ts(ta2,ta4,ta3) = 2delta8(tal,ta2) . 

ta3,ta4 

> the usual structure functions fs(tal,ta2,ta3) and ds(tal ,ta2,ta3) where all 
three indices vary from 1 to — 1. They are defined by the Lie algebra relation 

^Ts(tal,ta3,ta5) Ts(ta2,ta5,ta4) 

ta5 

2 

= — delta8(tal,ta2) delta3(ta3,ta4) 
Nc 

+ ^ ds (tal , ta2 , ta5) Ts (ta5 , ta3 , ta4) 

ta5 

+i^fs(tal,ta2,ta5) Ts (ta5 ,ta3 ,ta4) , 

ta5 

where the f s (tal ,ta2 ,ta5) are completely antisymmetric and the ds(tal,ta2, 
ta5) are completely symmetric. 

4.2.10 Output 

The output of the diagrams is decomposed in one or several form factors CI, C2, ... to 
reflect the Dirac 7 or Lorentz structure of the result. Therefore a list of form factors is 
returned. The last entry in this list is the defining equation for the form factors. The result 
is presented in the Maple text window of the main window and has the following structure: 

Gl := [CI = ...,..., Cn =..., CI ql(nul) q2(nu2) H h Cn G(nul,nu2) . . .1 

S v ' > „ ' 

the form factors Ci defining equation for form factors 

Let us consider for instance the self-energy of a vector boson, tfoops finds as Lorentz 
structure CI G(nul ,nu2) + C2 ql(nul) ql(nu2) and the output roughly looks like 

[CI = terml , C2 = term2 , CI G(nul ,nu2) + C2 ql(nul) ql(nu2)], 

where terml and term2 denote the form factors. In the case of a fermionic self-energy 
the result carries no Lorentz index. Instead, it is decomposed with respect to &*(1,0NE), 
&*(l,Dg5), &*(l,DgO) and &* (1 ,Dg5 ,DgO) and in the general case the result is of the 
form 

[CI = terml , C2 = term2 , C3 = termS , C4 = term4 , 
CI (1 &* DNE)+C2 (1 &* Dg5)+C3 (1 &* DgO) + C4 &*(1, Dg5, DgO)] . 

Here all 7-matrices are contracted. In vertex functions with a fermion line and a vector 
boson there is also a Dg(nul) appearing. In appendix [A] several examples are given in 
more detail. In the case of a scalar self-energy there appears only one form factor and the 
result reads 

[CI = terml , Cl] . 

The form factors can be written in several modes of evaluation which correspond to the 
different buttons of the diagram windows: 
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t> With ( Evaluate | or if no qualifier is added to the EvalGraph procedures the output - 
for simplicity - is returned in terms of unevaluated Oneloop and Twoloop functions 
which correspond directly to the OneLoop and TwoLoop functions explained in sect. 
O and O 



t> To get complete answers, ( Eval. Full] or in the case of EvalGraphl and EvalGraph2 the 
qualifier full forces 7(Coops to evaluate directly the Oneloop and Twoloop integrals. 
To be precise, you have to type 

EvalGraphm(n , (list) , full) ; 

where m and n are numbers. Now each form factor CI, . . . is a Laurent expansion in 
terms of the ultraviolet regulator e = (4 — D)/2, where the significant coefficients of 
this expansion - 0(e~ l ), O(e ) in the one-loop case, 0(e~ 2 ), 0(e^ 1 ), O(e ) in the 
two- loop case - are denoted in form of a list. The bosonic self-energy we introduced 
before now has the structure 

[CI = [e^ 1 -term, e° -term] , 
C2 = [e^ 1 -term, e° -term] , CI G(nul ,nu2) + C2 ql(nul) ql(nu2)] 

for one-loop contributions and 

[CI = [e~ 2 -term, e~ 1 -term, e° -term] , 
C2 = [e~ 2 -term, e~ l -term, -term] , CI G(nul ,nu2) + C2 ql(nul) ql(nu2)] 

for two-loop diagrams. e~ 2 -term and s~ 1 -term are meant to represent the divergent 
part (the coefficients of \ and -) whereas e°-term is describing the finite part of the 
corresponding form factor. For those two-loop topologies where no analytic result is 
known, the finite part contributes two entries to the output list. The first of them 
is a list itself which denotes an integral representation, the second contains the part 
which is analytically calculable: 

[CI = [e~ 2 -term, e~ 1 -term, numeric, e° -term] , 
C2 = [e~ 2 -term, e~ l -term, numeric, e° '-term] , 
CI G(nul,nu2) + C2 ql(nul) ql(nu2)]. 

A more detailed description of the integral representation - which was called here 



numeric - is given in sect. [4.4.3 



D> The integral representation is evaluated directly with the help of VEGAS if 
I Eval. Numeric | is used or if the qualifier full replaced by two lists which contain the 



input information for VEGAS: 
EvalGraphm(n, (list) , [mi,rri2] , {.n\,ri2\) \ 

Within an ordinary Maple session one has in addition to assign the variable NumPath 
as described in sect. [Oj. The variables m\,m2, n\, ni represent the grids where the in- 
tegrand is evaluated: first m\ iterations with m,2 points - itmx= m\ and ncall= rri2 
are passed to VEGAS - and then n\ iterations with n<i points - this corresponds to 
itmx= n\ and ncall= ri2- In the Xwindow interface these values can be changed 



with the entry Numeric of the Options menu (cf. sect. 4.6). Real and imaginary part 
are integrated one after the other. 
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The output has then the following structure: 

[CI = [e~ 2 -term, e~ 1 -term, [factor, [result, error]], e° -term] , 
C2 = [e~ 2 -term, e~~ l -term, [factor, [result, error]], e° -term] , 
CI G(nul,nu2) + C2 ql(nul) ql(nu2)] 

where factor denotes an analytically evaluated pre-factor of the integrand. The re- 
sult of the integral itself is called result here, error describes the uncertainty of the 
numerical integration which is returned by VEGAS. To be specific: the entire O(e ) 
contribution is obtained by 

factor x result + e°-term, 

the complete uncertainty is 

factor x error, 
where the uncertainty of e°-term is neglected. 

If not all parameters of the integrand carry numerical values the numerical integra- 
tion cannot be performed. Therefore the same output as for the qualifier full is 
returned. 

O With [ Evai. More ) or in the case of EvalGraphl, with the qualifier more instead of 
full, j(Coops evaluates even the 0{e) of one-loop integrals. This becomes relevant 
if the one-loop integral is part of a two-loop calculation, because then it may be 
multiplied with a divergent one- loop Z- factor. This changes the output of a bosonic 
self-energy: 

[CI = [e~ l -term, -term, e l -term] , 
C2 = [e~ 1 -term, e°-term, e 1 -term] , CI G(nul ,nu2) + C2 ql(nul) ql(nu2)]. 



As far as the masses and momenta which enter the EvalGraph functions are symbols the 
output is given algebraically, whereas if they are all assigned with numbers the result will 
be numerical. 

Generally, in the numerical case, the program sets the value of the imaginary part of 
the propagators, five digits higher than the numerical accuracy of the whole calculation, 
for example to 10~ 15 if one calculates with 20 digits. When reading the program Maple 
sets the number of digits to 20 - except when the user specified a higher accuracy before, 
for instance: 



Digits : =40; 



In the algebraic case just rho is returned. The output is facilitated by the usage of several 
abbreviations (4.22, 4.24, 4.25, 4.26| ) for certain one-loop contributions which are intro- 
duced in detail in the context of the OneLoop procedure in sect. |4.3.5 . The abbreviations 
are substituted if the evalRex command is used (cf. sect. 4.3.3j). 



4.2.11 Special relations 

If external momenta are zero or parallel, \Coops omits those Lorentz vectors which are 
no longer independent and therefore redundant. The number of form factors - and the 
dimension of the parallel space - may decrease. If this is the case the form factors will no 
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longer be expressed in terms of proper OneloopnPt and TwoloopnPtm functions, but with 
the help of its variants Oneloopfc_nPt and Twoloopfc_nPtm which are reduced in dimension 
of parallel space. The Oneloopfc_nPt functions correspond to the OneLoop&_nPt functions 
from sect. |4.3| . 

If for instance in the case of the two-point function the external momentum ql shall 
vanish, one has to make the following assignment before starting the calculation: 



qlO:=0; 

Now the bosonic self-energy reduces to 



[CI = terml, CI G(nul,nu2)] 



which means that one form factor vanished. The only remaining form factor is expressed by 
a Dneloopl_2Pt function - this is a two-point function with a parallel space of dimension 
zero. The notation "l_2Pt" reflects the fact that this function carries properties of the 
two-point function as well as of the one-point function. 

If in the case of a three-point function both external momenta ql and q2 have a com- 
mon rest frame one has to make the assignment 

q21:=0; 

Now the momentum q2 will disappear from the equation which defines the form fac- 
tors since it is no longer independent of ql ("zero recoil limit"). In addition the proper 
0neloop3Pt will switch to 0neloop2_3Pt functions - a mixing of two- and three-point 
functions. 

It is important here to make the assignment first and then to evaluate the diagram, 
because otherwise number and shape of form factors may be wrong. Moreover, if special 
relations exist which fix the variables of the EvalGraph functions or connect them with 
each other - like for instance the on-shell condition which connects masses and momenta 
- it is important first to declare these relations and then evaluate the functions. 

Otherwise, if the functions are evaluated first analytically and then numbers or re- 
lations between variables are substituted in the result, Maple has to do more work than 
necessary. Moreover, it might happen in some special cases that an error occurs - some- 
thing like "division by 0". Usually this sort of divergence is artificial and not an error 
%Coops is responsible for - at the moment of calculation the parameters which cause the 
error message were general analytical expressions. The problem can be overcome in ei- 
ther expanding in the critical parameters before the insertion of values, or - and this is 
recommended here - in calling the diagram window or the EvalGraph routines with the 
final values for the arguments from the very beginning, because the program was taught 
to detect a lot of this critical points and to avoid such artificial infinities in the result. 



4.2.12 Fast evaluation 

If the entry Oneloop Library of the Options menu is selected or - within an ordinary Maple 
session - the variable LibPath is assigned with some path, the procedures assume that 
a library of OneLoop functions is stored in the assigned directory - as described in sect. 
1.5 . The necessary integrals are then read from this directory - which of course is much 
faster (see sect. B.3). Otherwise, if LibPath is not assigned, the procedures are forced to 
recalculate the necessary OneLoop functions at any occurrence. 
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4.3 The OneLoop procedures 

4.3.1 Input 

There exist nine basic functions for the calculation of one-loop integrals. The first three 
functions 

> DneLooplPt(p, m, t) 

D> 0neLoop2Pt(p ,P_i_, <?io, mi, m 2 , h, h) 

> 0neLoop3Pt (p , Pi , P± , qio , <?20, 921 , mi , m 2 , m 3 , ti , t 2 , £3) 

correspond to the notation in parallel and orthogonal space variables. The next three 
functions 

> 0neLoopl_2Pt(p, mi, 7712, ii, t 2 ) 

[> 0neLoopl_3Pt(p, mi, m2, ffi3, ii, *2, £3) 

> 0neLoop2_3Pt(p ,P±, <?io, Q20, m 1 ,m 2 , m 3 ,ti,t2,t 3 ) 

are used in the case of reduced parallel space. All abovementioned functions require the 
tensor degree, the external momenta q n and the masses m n as input. The powers t n of 
the propagator terms in the denominator are optional. The remaining three functions 

D> OneLoopTenslPt(z, m, t) 

> 0neLoopTens2Pt(i, q\, mi,m 2 , ti,t 2 ) 

> 0neLoopTens3Pt (i , qi , q 2 , (92 — <?i ) 2 , mi , m 2 , m 3 , ti , t 2 , t 3 ) 

require the squared external momenta as input variables and return the complete rank-i 
tensor. The powers t n of the propagator terms in the denominator again are optional. 

4.3.2 Notation 

In detail the basic functions directly correspond to the following integrals: 

> One-point function: 





(4.5) 



DneLooplPt(p, m) = DneLooplPt(j), m, 1) 



(4.6) 




(4.7) 



DneLoopTenslPt(i, m) = OneLoopTenslPt(i, m, 1) . 



(4.8) 
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t> Two-point function: 




OneLoop2Pt(>o,P±, 9io, m 1 ,m 2 ,t 1 ,t 2 ) = £ {pop±){tl * 2) (<7io, m 1 ,m 2 ) 

= f dDl (l±) P± 

J [(l + qi) 2 -ml + ig]^ [f-m^ + ig]^ 

0neLoop2Pt(p ,P±, 9io, ^1,^2) = 0neLoop2Pt(p ,P±, Q10, mi, m 2 , 1, 1) 
0neLoopl_2Pt(p, mi,m2, ti,t 2 ) 

d 2 )5 



[Z 2 - mf + i^]* 1 [/ 2 - m 2 , + i£»]*2 
DneLoopl_2Pt(p, mi, 1712) = 0neLoopl_2Pt(p, mi,m2, 1, 1) 

DneLoopTens2Pt(i,qi,mi,m 2 ,ti,t 2 ) = B^: t . 2 J i (ql,mi,m 2 ) 



[(' + Qi) 2 - "if + ^l* 1 [I 2 — m| + i£>]* 2 
0neLoopTens2Pt(i, q\,mi,m 2 ) = 0neLoopTens2Pt(i, q 2 ,mi, m 2 , 1, 1) . 



Abbreviations: 



910 



/„2. j _ l -qi 
qi, to - 



z 2 -z 2 



(4.9) 
(4.10) 

(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 



t> Three-point function: 




0neLoop3Pt (p , Pi , RL , 9io , 920 , 921 , mi , m 2 , m 3 , h , t 2 , t 3 ) 

= cto>PW±Wut2,t 3 )( qi0j g20; q21 ,m u m 2 ,m 3 ) 



I 



d D l 



(lo) Po (h) Pl (l±) 



p± 



[(I + qi) 2 -m\ + igf 1 [(I + q 2 ) 2 -m\ + ig] t2 [I 2 - m§ + ig] t3 



(4.16) 
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0neLoop3Pt (p , p 1 , p ± , q 10 , q 20 , q 2 \ , m 1 , m 2 , m 3 ) 

= 0neLoop3Pt(p ,Pi,P±, 910,920, 921, mi, m 2 , m 3 , 1, 1, 1) (4.17) 
0neLoop2_3Pt(p ,£>±, 9io, 920, m 1 ,m 2 ,m 3 ,ti,t 2 ,t 3 ) 



[(! + 9i) 2 -m\ + igf 1 [(I + g 2 ) 2 - m\ + ig] t2 [I 2 - m§ + ig] t3 
0neLoop2_3Pt(p ,RL, 9io, 920, m 1 ,m 2 ,m 3 ) 

= 0neLoop2_3Pt(p ,P±, 9io, 920,mi,m 2 ,m 3 , 1,1,1) (4.18) 
0neLoopl_3Pt(p, mi,m 2 ,m 3 ,ti,t 2 , t 3 ) 

_ />, £>! 

7 [Z 2 - mf + i^]' 1 [Z 2 - m 2 + [J 2 - m| + ig)^ 

0neLoopl_3Pt(p, mi,m 2 ,m 3 ) = 0neLoopl_3Pt(p, mi,m 2 ,m 3 , 1, 1, 1) (4-19) 
0neLoopTens3Pt(i,qi,Q2, (92 - qi) 2 ,m 1 ,m 2 ,m 3 ,t 1 ,t 2 ,t 3 ) 

= C^ 3 \qlql(q 2 - Ql ) 2 , m u m 2 , m 3 ) (4.20) 

/"^D/ ' ' ' 

y [(Z + gi) 2 - mf + ig}^ [(I + q 2 ) 2 - m| + i£>]* 2 [Z 2 - m§ + i£>]* 3 
0neLoopTens3Pt(i, g 2 , g§, (92 — qi) 2 ,mi,m 2 ,m 3 ) 

= 0neLoopTens3Pt(i, g 2 , g|, (q 2 — gi) 2 , mi, m 2 , m 3 , 1, 1, 1) . (4.21) 
Abbreviations - using the auxiliary property (|4.2|): 



,2 91 * 92 U<?1 • 92) 2 2 

910 — V Ql ; 920 = — ^T" ! 921 = a/ 2 92 ; 



9i 



Z • gi , Z • g 2 



91 V-92 



, <<_l — y i o n 



z ± = wz 2 -z 2 -z 2 



(4.22) 



Our notation for the OneLoopnPt functions distinguishes between parallel and orthog- 
onal space which is reflected by the definitions of ( [4.15 ) and (4.22) for the momentum 



components. Please keep in mind the fact that in our notation the indices po,pi, ■ ■ ■ rep- 
resent the powers of the different components of the loop momentum Z which should not 
be mixed with Lorentz indices fii,u, 2 , ... of the corresponding tensor OneLoopTensnPt. 
If the powers t n of the propagator terms in the denominator are omitted like in ( [L6| , 



4.5, 4.10, 4.14, 4.17, 4.21) the program assumes the value 1 for each t n . 



4.3.3 Output of OneLoopnPt 

The results of the OneLoopnPt functions are Laurent expansions in terms of the ultraviolet 
regulator e. Therefore the output of the OneLoopnPt procedures consists of a list where 
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the significant coefficients - 0(e~ 1 ), 0(e ) - of this expansion are denoted. In two-loop 
calculations the term of 0(e 1 ) of the one-loop functions is also of interest. For that purpose 
this coefficient is also calculated if the qualifier more is used in the function call: 



Eingabe 


Ausgabe 


DtipT oonl P1~ I n 771 f) 


\?~~~ ^- fpVYYl - fp T77t\ 
C t/C'l III. C vC-l 1 1 v\ 


OneLooplPt(p, m, i,more) 


[e~ 1 -term, e^-term, e 1 -term] 


0neLoopl_2Pt(p, mi, m2, t±,t 2 ) 


[e~ 1 -term, e^-term] 


0neLoopl_2Pt(p, mi, m 2 , t±, t2, more) 


[e~ 1 -term, e^-term, e 1 -term] 


0neLoop2Pt (p , p± , qio , m 1 , m 2 , h , t 2 ) 


\£~ l -term, e^-term] 


0neLoop2Pt(p ,P±, qio,mi,m 2 , ti,t 2 ,more) 


[e~ 1 -term, e°-term, e 1 -term] 


0neLoopl_3Pt(p, mi, m 2 , m 3 , t±, t 2 , t 3 ) 


[e~ 1 -term, e°-term] 


0neLoopl_3Pt(p, mi, m 2 ,m 3 ,ti,t 2 , t 3 , more) 


\e~ l -term, e^-term, e 1 -term] 


0neLoop2_3Pt(p ,P±, <7io, 920, m 1 ,m 2 ,m 3 ,t 1 ,t 2l t 3 ) 


[e~ 1 -term, e°-term] 


OneLoop2_3Pt(po;P± 5 Qw, Q20, m i, m 2, r n-3,ti,t 2 , t 3 , more) 


[e~ 1 -term, e^-term, e 1 -term] 


0neLoop3Pt (p , p\,p±, qio ,q20,Q2i,mi,m 2 ,m 3 ,ti,t 2 ,t 3 ) 


[e~ 1 -term, e^-term] 


0neLoop3Pt (p , p 1 , p ± , q w , q 20 ,q 21 , 

mi,m 2 ,m 3 ,t 1 ,t 2 ,t 3 ,more) 


[e~ l -term, e°-term, e 1 -term] 



In the output e 1 -term is meant to represent the divergent part (the coefficient of ^) 
whereas e°-term is describing the finite part and e l -term the 0{e) contribution. 

As long as the arguments of the OneLoop functions are symbols the output is given 
algebraically, whereas of course numbers inserted for the arguments imply a numerical 
result. Like in the case of the EvalGraph routines, in numerical evaluations the program 
sets the value of g, the imaginary part of the propagators, five digits higher than the 
numerical accuracy of the whole calculation, for example to 10~ 15 if one calculates with 
20 digits. In the algebraic case just rho is returned. To get a compact result also in algebraic 
calculations, the output is written by using several abbreviations: 



R2exl(x, y) 



l- X - 



ln 1 



l- X - 



ln ( 1 + Jl - - ) +iir 

y 



(4.23) 



ln( — x) — ITT 



R2ex2(x,y) = 




+ 1 



In 1 + Jl 



(4.24) 



+ - (lny) 2 + 2 (lnx) 2 - 2 In a; In; 
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+ (In y — 2 In x 



1 + 4/1 



In 1 



y 



l- X - In l + 4 



ivr ^ X y [2 In 2 + ln(a - y)] 



-y 



R3ex2(x,y, z) = 2 In ( 1 

+2Li 2 (1 



r?(x,z) + 2 In (1 



+ 2Li 2 1 



R3ex3(x,y, z) = 4Si2 (l 



4Li 3 1 



+ 2 (lnz) 2 



+ 4Si2 1 



(4.25) 



4 Li, 1 



4r? he, - W Li 2 - + In 1 - - In - - + - 



1 



4r/h/, - MLi 2 -+lnl--ln--+- 



y 



y 



1 



In (1 
Infl- 



V 

z 

X 

z 

V 

z 



-4 lnz 
- 4 lnz 

+2Li 3 



Li 2 ( 1 - - ) + Li 2 ( 1 - - 



rj[x, ~)ln(l--)+7/(y, - I In ( 1 - ~ 



+ 



In 



z - y 

z — X 

g - y 

z — X 
12 



2 Li 



z 

(z - y)x 



In 



+2 (lnx — lny) 



g - y 

z — x 
Lio 



(z - x)y 
y — x z 
z ' z — X 

2 

-In 



+ 2Li 3 



(4.26) 



+ 2rj I y — x 



z — y \ I z — y z 
V 



(z - y)x 
(z - x)y 



z — X 



Lio - 



z z — X 



lull - -) (lnx- lny) 2 - 2C(3) - - (lnz) 



+2 



In x — In y In 

y 2 



(g ~ 
(z - x)y 



In 



(z — y)x\ ( y — x z 
(z - x)yy 77 



y z — x 



These functions are related to the corresponding 1Z functions^ [|T8[|T^]. They represent the 
coefficients of the Taylor expansion in e making use of (cf. |20|1 ) 



00 1 

cw = £ = 



k=l 



Li 2 (z) 
Li 3 (z) 



ln(l - s) 



ds 



Uo(s) 



ds 



The function R3ex3(x, y, z) assumes that x and y have an imaginary part of different sign which is 
always the case. 
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H21-I - r / — — — ds 



s 


17(0,6) = 2m[6(-Ima)e(-lmb)e(1m(ab)) - 6(lma)9(lmb)9(-lm(ab))] . 
The function 
evalRex(( expression)) ; 

substitutes the Rnexrn functions in (expression) by the corresponding Poly logarithms. 
4.3.4 Output of OneLoopTensnPt 

The functions OneLoopTensnPt return a full tensor as output. 

The output format looks similar to that of the EvalGraph routines. The following 
types of arguments - here only demonstrated for the two-point case - are allowed for all 
OneLoopTensnPt functions: 

> 0neLoopTens2Pt(i): 

The rank i tensor decomposition of the two-point function is given. The procedure 
returns a list consisting in the different coefficients, which are expressed in terms of 
the OneLoopnPt tensor integrals, and the defining equation for the coefficients, for 
instance in the case i = 2: 

0neLoop2Pt(0, 2) 0neLoop2Pt(0, 2) 

C2 1 — ~~ ~~ j C20 



3 - 2e ' q 2 (-3 + 2e) 

-f 0neLoop2Pt(2,0) , C20 ql(nul) ql(nu2)+C21 G(nul,nu2) 
T 

t> 0neLoopTens2Pt(i,full): 

The function inserts the results of the OneLoopnPt tensor integrals explicitly. Here 
again the case i = 2: 

[C21 = [e^ 1 -term, e° -term] , C20 = [e^ 1 -term, e° -term] , 

C20 ql(nul) ql(nu2)+C21 G(nul,nu2)] . 

t> 0neLoopTens2Pt(i, q{, m\,m2,t\,t<i): 

The function returns the same as 0neLoopTens2Pt(i, f ull), but expressed in the 
user-defined terms for q\,m\ and mi- Of course numerical values are also allowed. 
t\ and t<i are optional. 

t> 0neLoopTens2Pt(i, q{, mi,m2,ti,t2,^oxe): 

The function returns also the 0(e) contribution. Again, t\ and t 2 are optional: 

[C21 = [e~ l -term, e°-term, e 1 -term] , C20 = [e~ l -term, -term, e l -term] , 
C20 ql(nul) ql(nu2)+C21 G(nul,nu2)] . 



The C20, C21, 



correspond to the coefficients of the Passarino-Veltman procedure. 
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4.3.5 Infrared divergences 

The program is regulating the infrared divergence of the three-point function as well - 
which may occur for instance if the momenta are on-shell and one mass is set to 0. There is 
no distinction made between UV and IR divergences - that means that there is no infrared 
dimension parameter ErR. Both kinds of divergences are described by e = (4 — D)/2. 



4.3.6 Special relations 

If special relations exist which fix the variables of the OneLoop functions or connect them 
with each other - like for instance the on-shell condition - it is important first to declare 
these relations and then to call the functions - as already stressed in the case of the 
EvalGraph routines. 

If the functions are evaluated first analytically and then numbers or relations between 
variables are substituted in the result, it might happen in some special cases that an error 
occurs - something like "division by 0". Usually this sort of divergences is artifical due 
to the way the analytical result was written before the insertion (cf. sect. 4.2.11] ). 



4.4 The TwoLoop procedures 
4.4.1 Input 

Unlike for the one-loop functions on the two-loop level for each n-point case exist different 
types of integrals - which correspond to different topologies with the same number of 
external legs. To distinguish between the different n-point functions with the same n the 
TwoLoop procedures acquire an additional identification number: 

> TwoLoop2Ptl(p ,P±, r , r±, s, q w , mi, m 2 , m 3 , m<t, m 5 , h, t 2 , t 3 , £ 4 , t 5 ) 

> TwoLoop2Pt8(p ,P±, r , r±, s, qw, mi, m 2 , h, t 2 ) 

t> TwoLoopTens2Ptl(ii, i2,qi,mi,m2, m 3, m i, m 5,ti,t2,t3,t4:, t^) 



> TwoLoopTens2Pt8(ii, i2,<li, m i, m 2,ti,t2)- 

Like in the one-loop case each function needs the tensor degree, the external momenta q n 
and the masses m n as input. The powers t n are again optional. 

The functions TwoLoopTensnPtm require the squared external momenta as input vari- 
ables and return the complete tensor, whereas the functions OneLoopnPtm correspond to 
the notation in parallel and orthogonal space variables. 
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4 Detailed description 



4.4.2 Notation 

The functions correspond to the following integrals: 
t> Two-point functions: 

Topology 1: 




l + qi 



TwoLoop2Ptl(p ,RL, r , rj_, s, q w , mi, m 2 , m 3 , m 4 , m 5 , h,t 2 ,t 3 , t 4 , t 5 ) 



J d D i J d D k 



[(l + Qi) 2 ~mj + ig]^ [I 2 -m\ + ig] t2 [l 2 - m\ + i£>]*3 

1 



[(i + fc) 2 -m\ + ig] u [k 2 - m\ + i^ps 
TwoLoopTens2Ptl(«i, Z2, qi,mi,m 2 ,m^, m.4, ms, ti, *2, *4, is) 
= fd°l j d D k 



[(' + <?i) 2 ~m\ + igY 1 [I 2 - to 2 , + ig]^ [I 2 - m§ + zp] 



[(i + fc) 2 -m\ + ig] u [k 2 - m\ + ig}^ ' 



*3 



Topology 2: 




TwoLoop2Pt2(p ,_p_L, r , r±, s, q 10 , m 1 , m 2 , m 3 , m 4 , m 5 , h,t 2 ,t 3 , i 4 , t 5 ) 

(Z ) w (fc ) ro (fc±) r - L 



= fd°lj d D k 



(4.27) 



(4.28) 



(4.29) 



[(' + <?i) 2 -m\ + ig]^ [I 2 - m 2 , + ig] t2 [(I + k) 2 - m\ + 



*3 



[(fc - <?i) 2 -m\ + ig] u [k 2 -m\ + ig\ t5 
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TwoLoopTens2Pt2(i 1 , i 2 , Qi, mi, m 2 ,mz, 1714,1715, t\,t 2 , £3, i 4 , t 5 ) 

= !d D il d D k ^ " ' ^ " ' kui2 



(4.30) 



[(I + gi) 2 - m 2 + [Z 2 - m 2 + ig]ti [(I + fc) 2 - m\ + i e ]*s 

1 



[(fc - qi) 2 - m| + i£>]* 4 [fc 2 - raf, + i£>]*5 



Topology 3: 




TwoLoop2Pt3(p ,RL, r , r±, s, qio,m 1 ,m 2 , m 3 , m 4 , h,t 2 ,t 3 , t 4 ) 



(/o) po (fc ) ru 
[(i + Qi) 2 -m\ + igY 1 [I 2 - m| + i^]* 2 

1 



[(Z + fc) 2 - m§ + i£»]* 3 [fc 2 -m\ + ig] u 
TwoLoopTens2Pt3(ii, 9i> mi,m 2 , 1713,1714, t±, t 2 , t 3 , i 4 ) 

lfll' " Vii ^1 ' ' ' k"» 2 
[(* + <?i) 2 -m\+ ig]^ [I 2 - m| + i^]* 2 

1 



/ rfD '/ 



d°A; 



[(Z + fe) 2 - mjj + i£»]* 3 [fc 2 - m 2 + i£»]' 4 



Topology 4: 





« + fe 






rm. 




> 4 


m 3 


-* > 




> 

, k 

^v. mi 

l + qi 





(4.31) 



(4.32) 
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4 Detailed description 



TwoLoop2Pt4(p ,P_L, r , r±, s, q w ,m 1 ,m 2 , m 3 , h,t 2 ,t 3 ) 

(/ )«> {i ± y>± (k y° (k ± y± ( z y 



(4.33) 



= Jd D l J d D k 



[(I + Qi) 2 -m\ + igY 1 [(I + k) 2 - m 2 , + ig] t2 [k 2 - m§ + ig] t3 



TwoLoopTens2Pt4(zi, Z2, q\, mi, m 2 , m 3 , ti,t 2 ,ts) 



(4.34) 



/*"'/' 



[(* + <?i) 2 -m\ + ig]^ [(I + k) 2 - m?> + ig] t2 [k 2 - m\ + ig]^ 



Topology 5: 



Qi 




<7i 



l + qi 

TwoLoop2Pt5(p ,RL, r ,rj_, s, q w , mi, m 2 , m 3 , m 4 , ti, t 2 , *3, *4) 



(jo)w (Z ± )^ (fcoT° (fcjj^ 
[(Z + Qi) 2 -m\ + ig]*! [I 2 -m% + ig]^ 

1 



[Z 2 - m§ + i£>]* 3 [A; 2 -m\ + ig] 1 * 
TwoLoopTens2Pt5(ii, i 2 , q\, mi, m 2 , m 3 , 7714, ti,t 2 ,t 3 , £4) 

[(Z + <7i) 2 - mf + s^]* 1 [Z 2 - m 2 , + i^]* 2 

1 



[Z 2 - m§ + i£»]* 3 [fc -m\ + ig] u 



Topology 6: 




I + 91 



k — qi 



(4.35) 



(4.36) 



4 Detailed description 



55 



TwoLoop2Pt6(p ,RL, r , rj_, s, gi , mi, m 2 , m 3 , m 4 , <i, i 2 , t 3 , 



= d D l d D k 



(Zq)PO (/ ± )px (fc ) r ° (fcjj^ 
[(I + Qi) 2 -m\ + ig] h [I 2 -m\ + ig] t2 

1 



[(k — qi) 2 - m§ + i£»]' 3 [fc -m\ + z£>] 
TwoLoopTens2Pt6(«i, ^2, 9i, mi, 1712,1713,1714, ti, <2, £3, £4) 



[(* + <?i) 2 -m\+ igY 1 [I 2 -m\+ igY 2 



(4.37) 



(4.38) 



Topology 7: 



<7i 



[(fc — <?i) 2 - m§ + i£>]' 3 [fc 2 -m\ + igY* 




<7i 



l + qi 



TwoLoop2Pt7(p ,P±, r , r±, s, q w ,m 1 ,m 2 , m 3 , ti, t 2 , t 3 ) 



d D l d D k 



(io) po (i±) p± (ko) ro (k ± y± ( z y 



1*1 r/ 2 _ ml 



ISn, 



[(I + Qi) 2 ~ + igY 1 [I 2 - m 2 , + ig\ t2 [k 2 - m§ + igY 3 



(4.39) 



[(I + Qi) 2 ~ml + igY 1 [I 2 - m\ + igY 2 [k 2 - m\ + ^]*» 
TwoLoopTens2Pt7(ii, i 2 , q\, m\,m 2 ,m 3 ,t\,t 2 , i 3 ) (4.40) 
= JdPlJ d D k 
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4 Detailed description 



Topology 8: 




l + qi 

TwoLoop2Pt8(j> ,RL, r , r±, s, q w , mi, ma, h, i 2 ) (4.41) 
d D l (d D k {lo)P ° {l±)P± {koY ° {k±Y± {Z)S 



[{I + Qi) 2 -m\ + igY 1 [k 2 - m\ + ig] t2 



TwoLoopTens2Pt8(ii, i2, <?i, mi, m2, ii, £2) (4.42) 



d D l ld D k 



[(I + qi) 2 -mf + iqY 1 [k 2 -m% + ig] t2 ' 

Like in the one-loop case for each function holds 

TwoLoop2Ptn(po,P±,r ,r_L,s,gi ,mi,. . .,m k ) 

= TwoLoop2Ptn(p ,P±, r , r±, s, q w , mi, ... , m k , 1, . . . , 1) (4.43) 
TwoLoopTens2Ptn(ii, 72, Qi, tni, ■ ■ ■ , f^k) 

= TwoLoopTens2Ptn(7i, ?2, Qi, m i) • • • > m fcj 1, • • ■ , 1) - (4.44) 

This means that if the powers t n of the propagator terms in the denominator are omitted 
the program assumes the value 1 for each t n . 



The abbreviations ( 4.15 ) we made in the one-loop case for parallel and orthogonal 
space notation are still valid. 

There are several relations between different two-point functions. Some of them we 
would like to mention here: 

TwoLoop2Pt 1 (p , p ± , r ,r ± , s , q i0 , mi , m 2 , m 3 , m 4 , m 5 ) 

: 2 [^oLoop2Pt3(p , p ± ,r ,r ± ,s,q 1 o,m 1 ,m 3 ,m 4 ,,m 5 ) (4.45) 



ms — m 



2 



-TwoLoop2Pt3(p ,P±, r 0) r ± , s, q 10 , m x , m 2 , m 4 , m 5 )] 



TwoLoop2Pt5(p ,P_L, r , r_L, s, <?i , mi, m 2 , m 3 , m 4 ) 

= —9 o [TwoLoop2Pt7(p ,P±,ro,r_ L ,s,gio,mi,m 3 ,m4) (4.46) 

m 3 ~~ m 2 

-TwoLoop2Pt7(p ,p ± , r , r ± , s, g l0 , ra 1} m 2 , m 4 )] . 
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Both relations are only valid for m 2 7^ m 3 , otherwise one has 

TwoLoop2Ptl(p ,P±, r , r±, s, q 10 , mi, m 2 , m 2 , m 4 , m 5 , h,t 2 ,t 3 , t 4 , t 5 ) 

= TwoLoop2Pt3(p ,P±, r , rj_, s, gi , mi, m 2 , m 4 , m 5 , *i,t 2 + *3, *4, *s) (4.47) 

TwoLoop2Pt5(p ,RL, r , r_L, s, gio, m\,m 2 ,m 2 ,m<i, h,t 2 , t 3 , U) 

= TwoLoop2Pt7(p ,P±, r ,r_L,s, #10, mi, m 2 , m 4 , t l5 i 2 + <3,*4)- (4.48) 

The two-point functions numbered with 1-4 are nonfactorizable whereas the functions 5-8 
are products of one- loop integrals, for instance 

TwoLoop2Pt6(p ,RL, r , r±, s, q w , mi, m 2 , m 3 , m 4 , t x , t 2 , t 3 ,U) (4.49) 
0neLoop2Pt(p ,P±, Qio, m 1 ,m 2 ,t 1 ,t 2 ) 



-2' 2' 
D-l 



S.) 

y 2 > 2> 

x OneLoop2Pt(ro, r±, — qio, m 3 , m 4 , t 3 , t 4 ) if s even 
if s odd 



TwoLoop2Pt7(p ,RL, ro, r_L, s, qio,mi,m 2 ,m 3 ,t 1 ,t 2 , t 3 ) (4.50) 
(if) , ^ (D-l)---(D-3 + r ± ) r ! 



(^i,|) v ' D ... {D - 2 + r + r ± ) (a )!2 ? 

x 0neLoop2Pt(p ,P±, Qio, rni,m 2 ,ti,t 2 ) 

x OneLooplPt(ro + r_|_, 7713, £3) if s, ro, rj_ even 

else . 



4.4.3 Output of TwoLoopnPtm 

The output - like in the one-loop case - is a list which contains the relevant coefficients of 
the Laurent expansion in e. Ultraviolet divergences at the two-loop level occur as terms 
of 0(e~ 2 ) and 0(e~ 1 ). Therefore the notation of the output starts with two divergent 
coefficients. 

For the two-point functions 1-4 the finite part cannot be given completely analytically 
for all mass cases. Therefore these functions have two entries in the output list which 
describe the finite part. The first one is a list itself denoting an integral representation, 
the second one gives the analytically calculable terms. 

The integral representation is evaluated directly with the help of VEGAS if two argu- 
ments are added. Both arguments must be lists which contain each two numbers. These 
numbers carry the input information for VEGAS. Within an ordinary Maple session one 
has in addition to assign the variable NumPath as described in sect. [4.6| . 



input 


output 


TwoLoop2Ptm(p ,P±, ^0, r±,s, q w , 
mi, . . . ,m k ,ti, ... ,t k ) 


[e~' 2 -term, e~ 1 -term, 

numeric, e°-term] 


1^oLoop2Ptm(po,p ± ,r ,rj_,s,qio,mi, . . . ,m k , 
ti, ...,t k , [h,l 2 ], [ni,n 2 ]) 


[e~ 2 -term, e~ 1 -term, 

[factor, [result, error]], e° -term] 
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4 Detailed description 



numeric describes the integral representation. This list consists of the following entries 

[factor, integrand, x=- infinity . .infinity, y=-infinity. .infinity] 

where integrand is the integrand of the integral representation which has to be integrated 
over x and y in the interval (— 00,00). factor is an analytical factor which multiplies the 
integral representation. 

The variables h,h, n i, n 2 represent the grids where the integrand is evaluated: first 
l\ iterations with I2 points - itmx= l\ and ncall= l 2 are passed to VEGAS - and then 
n\ iterations with n 2 points - this corresponds to itmx= n\ and ncall= n^. Real and 
imaginary part are integrated one after the other. 

The other entries of the output list have the same behaviour as in the one-loop case: 
As long as the arguments of the TwoLoop functions are symbols the output is given al- 
gebraically, whereas numbers inserted for the arguments imply a numerical result. In the 
numerical case, the program sets the value of g, the imaginary part of the propagators, 
five digits higher than the numerical accuracy of the whole calculation, for example to 
10~ 15 if one calculates with 20 digits. In the algebraic case just "rho" is returned. 

In the case where all internal particles are massless the result is of course known 
analytically. In this case numeric is just [0] . 

The two-point functions 5-8 are products of one-loop integrals and therefore evaluated 
completely analytically. The output list has no entry for numerical evaluation and looks 
simply like: 



input 


output 


TwoLoop2Ptm(p ,^±, r , r±,s, q w , 
mi, . . . ,m k ,h, . . . ,t k ) 


[e~ 2 -term, e~ l -term, e°-term] 


TwoLoop2Ptm(po,P±,n),r_L,s,q<io,mi, • • • ,m k , 
ti, ...,t k , [h,h], [ni,n 2 ]) 


[e~ 2 -term, e~ l -term, e°-term] 



The optional input of two lists for the VEGAS input data doesn't affect the output since no 
numerical integration is performed here. The 0{e 1 ) is not relevant in two-loop calculations. 
Therefore the qualifier more does not exist for two-loop functions. 

The divergent part of two-loop integrals is closely related to one-loop integrals. There- 
fore the same abbreviations ( 4.23 , 4.24f) as for one-loop integrals are appearing. They are 
substituted if the evalRex command is used (cf. sect. 4.3.3). 



4.4.4 Output of TwoLoopTensnPtm 

The functions TwoLoopTensnPtm return a full tensor as output. The same types of ar- 
guments are allowed as for the corresponding OneLoopTensnPt functions described in 



sect. 4.3.4 - except the qualifier more. 



4.5 The OneLoopLib procedures 

For different calculations of the same OneLoop function it is of course rather inconvenient 
to start the whole program again. Therefore our package includes the functions 

[> 0neLoopLib2Pt (i,j) 

> 0neLoopLib3Pt (i,j) 



4 Detailed description 
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which generate a library of all 0neLoop2Pt and 0neLoop3Pt functions respectively up to 
the tensor rank i. All powers of the denominators from 1 to j are considered. The second 
parameter j is optional. If it is omitted the procedures assume the value 1 for j, so that 
no powers of denominators higher than one are calculated. These functions are then read 
by the proper procedures described before, so that they in any case speed up. 

If you use the OneLoopLib functions within an ordinary Maple session and want to 
write this library in any other than the actual directory you have to assign the variable 
LibPath: 



LibPath : =(path) ; 

Now the library will be written to (path). The procedures will only look for the library if 
the entry Oneloop Library of the Options menu is selected or if LibPath is assigned. If the 
procedures search for an integral which is not contained in the library it will be calculated 



and stored in the library on the fly. Each integral corresponds to one file. In appendix B.2 
we list the convention for the file names. 



4.6 Numerical integration with VEGAS 

For numerical integration a parallelized implementation of the VEGAS algorithm for C++ 
invented by R. Kreckel is provided with tfoops. To use the routines within an ordinary 
Maple session one has to declare the variable NumPath: 

NumPath:=(pa£/i) ; 

(path) describes the directory where the integration routines are located. If Troops is in- 
stalled properly this is the subdirectory cxx of the ^loops directory. 

The routines are customized for those TwoLoop procedures which cannot be calculated 
analytically. The most elegant way to perform those numerical integrations is described 
in sect. 4.2. 10| and 4.4.3. In addition there is also the possibility to perform an integration 



separately. For that purpose there exists the command 
NumlntC ( (integrand) , [mi, 771,2] > [14,^2]) ! 

(integrand) denotes the expression which has to be integrated numerically, xloops expects 
the integration variables to be named x and y - as they are called by the EvalGraph2 and 
TwoLoop procedures. The variables mi,m2,ni,n2 represent the grids where the integrand 
is evaluated: first mi iterations with 777,2 points - itmx= mi and ncall= 7712 are passed 
to VEGAS - and then ni iterations with ri2 points - this corresponds to itmx= rii and 
ncall= 71-2- Real and imaginary part are integrated one after the other. If only the real 
part is of interest one has to specify 

NumlntC ( (integrand) , [mi, 771,2] , [771,712] ,Re) ; 

if only the imaginary part shall be evaluated one writes 

NumlntC ( (integrand) , [mi, m<i~\ , [771,772] ,1m) ; 



If NumlntC is called only with the first argument sloops takes the values mi = 20, 7772 = 
1000, m = 5 and n 2 = 10000. 
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4 Detailed description 



The output has the following format: 

[(result), (error)} . 

(result) describes the - possibly complex - numerical result, (error) is the uncertainty (la) 
of the numerical integration which is returned by VEGAS. 

NumlntC performs the necessary internal steps („ compile, link") and starts the C++ 
program as subprocess. Maple waits until the subprocess terminates and finally re-reads 
the result. 

The numerical integration with VEGAS should be used only once by each user. Other- 
wise there will occur conflicts because several numerical integrations would use the same 
C++ files. 



Vegas Options 



1st Integration 



itniK (max. number of it.) 20 



ncall (number of calls) |1 000 



2nd Integration 



itmx (max. number of it.) [s~ 



ncall (number of calls) 10000 



Ok 



Figure 4.7: Options for VEGAS 



The entry Numeric of the Options menu enables the user to change the parameters for 
VEGAS. These parameters are set in a window displayed in fig. 4/7. 



5 Conclusion and outlook 



Our aim was to introduce a program package that gives any computer user the possibility 
to calculate particle processes without much effort in manpower and hardware. At this 
stage of development the program is still understood to be under construction. In fact we 
are planning to incorporate several additional features: 

D> One-loop four-, five- and six-point functions. 

D> Two-loop three- and four-point functions. 

> Calculation of complete processes, generation of all diagrams which contribute to a 
process in a given order. 

> Incorporation of other models, for instance supersymmetric theories. 

> An interface for adding other models. 

> A plotting device for automatical plotting of results. 

> Automatical renormalization, determination of Z factors and renormalized functions. 
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Appendix 



A Examples 



A . 1 Installation 

A typical run of the configure script looks like the following: 

. /configure 

This is the XLOQPS installation script ! 

Which Release of Maple V are you using [1/2/3/4] ? 3 

Enter the name of the WWW Browser executable, 
which XLOOPS should use ! netscape 

Do you have Isolatin-1 character encoding [y/n] ? y 



If you feel unsure, answer the following question with 'n' . 

XLOOPS can work in single processor mode on multi processor machines as well. 
Do you have multiprocessing (and POSIX thread library installed) [y/n] ? y 
How many processors has your computer ? 2 

On which platform are you installing xloops ? 
Linux [1] Digital OSF [d] SUN Solaris [s] IBM AIX [a] Other Unix [o] 1 

Instead of creating the library you can also get a prebuilt library 
from the XLOOPS homepage and unpack it . 
Shall I create the library of one-loop functions 
(takes approx. 1 hour) [y/n] ? n 

Now testing for Maple output bug while reading the library ! 

IWI Maple V Release 3 (Universitaet Mainz) 

• _l\l l/l_. Copyright (c) 1981-1994 by Waterloo Maple Software and the 
\ MAPLE / University of Waterloo. All rights reserved. Maple and Maple V 

< > are registered trademarks of Waterloo Maple Software. 

I Type ? for help. 

> read( 'loops. ma' ) ; 
Warning: new definition for &* 



»»»»»»»»»»»»»»»»»» This is XLOOPS Version 1.0 ««««««««««««««««««« 

by Lars Briicher, Johannes Franzkowski, Dirk Kreimer 

©1993-1997 Johannes Gutenberg Universitat Mainz, Germany 

bytes used=1001220, alloc=851812 , time=0.03 
bytes used=2011736, alloc=1310480, time=0.03 
> LibPath := ' . ./lib/' ; 

LibPath := ../lib/ 
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A Examples 



> EvalGraphl (3 , [higgs , higgs , higgs , higgs , higgs , higgs] , full) ; 
bytes used=3011944, alloc=1834672 , time=0.03 

> quit ; 

bytes used=3662072, alloc=1834672 , time=0.03 

Did you see any result of the Maple command 'Eval Graph' used above [y/n] ? n 
System configuration completed, doing now some checks ! 
This may take a while ! 
Creating the index file for tcl ! 
Now testing the Makefile for numerical integrations generated for your System ! 
rm -f main f unct . o numint.o divO.o vegas.o gauss. o utils.o 

g++ -DGNU_COMPLEX -03 -f omit-f rame-pointer -f inline-functions -m486 -malign- 

functions=4 -malign-jumps=4 -malign-loops=4 -f expensive-optimizations -funroll-1 
oops -f fast-math -c f unct . cxx 

g++ -DGNU_C0MPLEX -03 -f omit-f rame-pointer -f inline-functions -m486 -malign- 

functions=4 -malign-jumps=4 -malign-loops=4 -f expensive-optimizations -funroll-1 
oops -f fast-math -c numint.cxx 

g++ -DGNU_C0MPLEX -03 -f omit-f rame-pointer -f inline-functions -m486 -malign- 

functions=4 -malign-jumps=4 -malign-loops=4 -f expensive-optimizations -funroll-1 
oops -f fast-math -c divO.cxx 

g++ -DGNU_C0MPLEX -03 -f omit-f rame-pointer -f inline-functions -m486 -malign-fun 
ctions=4 -malign-jumps=4 -malign-loops=4 -f expensive-optimizations -funroll-loop 
s -f fast-math -D_REENTRANT -DPVEGAS_P0SIX -c vegas . c 

g++ -DGNU_C0MPLEX -03 -f omit-f rame-pointer -f inline-functions -m486 -malign- 

functions=4 -malign-jumps=4 -malign-loops=4 -f expensive-optimizations -funroll-1 
oops -f fast-math -c gauss . cxx 

g++ -DGNU_C0MPLEX -03 -f omit-f rame-pointer -f inline-functions -m486 -malign- 

functions=4 -malign-jumps=4 -malign-loops=4 -f expensive-optimizations -funroll-1 
oops -f fast-math -c utils.cxx 

g++ -o main f unct . o numint.o divO.o vegas.o gauss. o utils.o -lm -lpthread mai 
n.cxx 

Division by Errors will be ignored 
Test: 1/0 

Initializing SR-sequences with seed 876920052 
Input parameters for vegas: ndim= 2 ncall= 968 2 thread(s) 

ittot= 1 itmx= 5 22" 1 hypercubes 

nprn= ALPH= 1.50 

mds= 1 nd= 200 npg=2 
xl[ 1]= xu[ 1]= 1 

xl[ 2]= xu[ 2]= 1 

iteration no. 1 : integral = -1 . 077957e-06 +/- 9.2e-07 

all iterations: integral = -1 . 077957e-06+/- 9.2e-07 chi**2/IT n = 
iteration no. 2 : integral = 1.703926e-07 +/- 7.1e-07 

all iterations: integral = -2 . 975584e-07+/- 5.6e-07 chi**2/IT n = 1.2 
iteration no. 3 : integral = 1.407306e-09 +/- 4.3e-07 

all iterations: integral = -1 . 079265e-07+/- 3.4e-07 chi**2/IT n = 0.67 
iteration no. 4 : integral = -6 . 310369e-08 +/- 4e-07 

all iterations: integral = -8 . 938296e-08+/- 2.6e-07 chi**2/IT n = 0.45 
iteration no. 5 : integral = -2 . 702396e-07 +/- 3.2e-07 

all iterations: integral = -1 . 618004e-07+/- 2e-07 chi**2/IT n = 0.38 
Input parameters for vegas: ndim= 2 ncall= 9800 2 thread(s) 

ittot= 1 itmx= 3 70~1 hypercubes 

nprn= ALPH= 1.50 

mds= 1 nd= 200 npg=2 
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: integral 
integral = 
: integral 
integral = 
: integral 
integral = 
Total Number of ignored Divisions by 0: 
Input parameters for vegas : ndim= 2 



iteration no. 
all iterations: 
iteration no. 
all iterations: 
iteration no. 
all iterations: 



xl[ 1]= xu[ 1]= 1 

xl[ 2]= xu[ 2]= 1 

-4.470697e-08 +/- 6.2e-08 
4.470697e-08+/- 6.2e-08 chi**2/IT n = 

2.028027e-08 +/- 5.3e-08 
6.763662e-09+/- 4e-08 chi**2/IT n = 0.63 

-2.663549e-08 +/- 5.1e-08 
1.435286e-08+/- 3.2e-08 chi**2/IT n = 0.36 


ncall= 968 



ittot 

nprn= 

mds= 1 
xl[ 1] 
xl[ 2] 



2 thread (s) 
22" 1 hypercubes 



: integral 
integral = 
: integral 
integral = 
: integral 
integral = 
: integral 
integral = 
: integral 
integral = 
Input parameters for vegas : 



iteration no. 
all iterations : 
iteration no. 
all iterations: 
iteration no. 
all iterations : 
iteration no. 
all iterations: 
iteration no. 
all iterations: 



iteration no. 
all iterations: 
iteration no. 
all iterations : 
iteration no. 
all iterations: 



1 itmx= 5 
ALPH= 1.50 

nd= 200 npg=2 

xu[ 1]= 1 
xu[ 2]= 1 
7.704579e-06 +/- 1.3e-06 
7.704579e-06+/- 1.3e-06 chi**2/IT n = 

8.286222e-06 +/- 6.8e-07 
8.160574e-06+/- 6.1e-07 chi**2/IT n = 0.16 
= 9.681104e-06 +/- 6.2e-07 
8.902334e-06+/- 4.3e-07 chi**2/IT n = 1.6 

8.987195e-06 +/- 5e-07 
8.938647e-06+/- 3.3e-07 chi**2/IT n = 1.1 
= 8.732733e-06 +/- 3.7e-07 

8.84922e-06+/- 2.5e-07 chi**2/IT n = 0.85 
ndim= 2 ncall= 9800 
ittot= 1 itmx= 3 
nprn= ALPH= 1.50 

mds= 1 nd= 200 npg=2 
xl[ 1]= xu[ 1]= 1 

xl[ 2]= xu[ 2]= 1 

: integral = 8.890832e-06 +/- le-07 

integral = 8 . 890832e-06+/- le-07 chi**2/IT n = 
: integral = 8.87138e-06 +/- 5.3e-08 

integral = 8 . 875575e-06+/- 4.7e-08 chi**2/IT n = 0.029 
: integral = 8.798005e-06 +/- 4.7e-08 

integral = 8 . 836898e-06+/- 3.4e-08 chi**2/IT n = 0.68 



2 thread (s) 
70" 1 hypercubes 



Total Number of ignored Divisions by 0: 



Make seems to work, so Makefile is correct ! 



If you had problems please edit the Makefile in the subdirectory ./cxx 
until it compiles correctly. Please report your system configuration and 
Makefile changes to xloops@thep.physik.uni-mainz.de ! 

Do you agree sending the XLQ0PS team a mail [Y/n] ? Y 

XL00PS installed successfully, invoke with ./xloops ! 
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A. 2 Examples with the Xwindows front end 

To demonstrate the input and output of tfoops some typical examples are given. The 
Feynman diagrams in this chapter are Postscript pictures produced with \Coops. 

A. 2.1 Output of one-loop integrals 
Analytic examples 



At first we want to show the different ways of evalu- 
ating the diagrams. As an example we have chosen a 
simple self-energy diagram. The advantage of this dia- 
gram is the relatively short result. The diagram on the 
right shows the process, a nondiagonal Z self-energy 
with a W bubble. 




With [ Evaluate | the following result is obtained: 
Gl := [ 

2 2 2 2 2 2 2 

CI = 1/16 I e cos(tw) (1 + eps Ln(4 Pi MU ) + 1/2 eps Ln(4 Pi MU ) ) 

/ 4 
(61-41 eps) OnelooplPt(0, Mw) / (sin(tw) Pi ), 

/ 

CI G(nul, mi2)] 

As expected, the structure of form factors, indicated by the last element of the list above, 
is just proportional to g^ u . The pre-factor CI is determined by an equation, which shows 
the corresponding Oneloop functions. In this example just OnelooplPt (0 , Mw) appears. 
As pre-factor the coupling e, the Weinberg angle tw and the renormalization parameter 
MU occur. 



With Eval. Full j the OneLoop function can be evaluated further: 
GF2 := [ 

2 2 
I e cos(tw) Mw 

CI = [- 3/8 , 

2 

sin(tw) Pi 

2 2 2 2 2 2 

1/16 I (- 6 Ln(4 Pi MU ) Pi Mw + 4 Pi Mw 

2 2 2 2 

- 6 (- Ln(Mw - I rho) - Ln(Pi) - gamma + 1) Pi Mw ) e cos(tw) 

/ 4 

/ (sin(tw) Pi ) ] , 

/ 



CI G(nul, nu2)] 
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Now the constant CI is evaluated to a list, which elements represent the coefficients of the 
Laurent series in the dimensional regulator e = (4 — D)/2. The first element of the list 
shows the divergence of the integral, the second element the convergent part. 

Additionally the next order in dimensional regularization can be evaluated. Clicking 
on ( Eval. More] shows the result: 

GM3 := [ 

2 2 222222 2 

I e cos(tw) Mw I (- 6 Ln(4 Pi MU ) Pi Mw + 4 Pi Mw - 6 7,1) e cos(tw) 

CI = [- 3/8 , 1/16 , 

2 4 

sin(tw) Pi sin(tw) Pi 

2 2 2 2 2 2 

1/16 I e cos(tw) (- 3 Ln(4 Pi MU ) Pi Mw 

2 2 2 2 2 2 

+ Ln(4 Pi MU ) (4 Pi Mw - 6 7.1) - 6 (1/2 Ln(Mw - I rho) 

2 2 

- (- Ln(Pi) - gamma + 1) Ln(Mw - I rho) + 1/2 Ln(Pi) 

2 2 2 2 

- (- gamma + 1) Ln(Pi) + 1/12 Pi + 1/2 gamma - gamma + 1) Pi Mw +4 7,1 

I 4 
) / (sin(tw) Pi ) 
/ 

]. 

CI G(mil, nu2)] 

2 2 2 

7,1 : = (- Ln(Mw - I rho) - Ln(Pi) - gamma + 1) Pi Mw 

In the above result the list for CI was just enhanced by an additinal element, the 0(e l ). 



Numerical examples 

In the following subsection we want to show some numerical examples, where different 
kinds of form factors occur. The values for couplings and particle masses where inserted 
with Insert Particle Properties from the Options menu. The momentum was set to qlO : =MzO ; 
and the renormalization parameter was chosen to be MU : =1 ; . All results are given in natural 
units H = c = 1 and the energy scale GeV. 




(qio.o) > (qi°,Q) > (qio,Q) (qio.Q) (qio,o) (qio.o) 

(a) (b) (c) 
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The results, when evaluated with ( Evai. Full ) , read: 
(a) 

GF4 := [ 

CI = [ - 2143.7993127296539430 I, 

-15 

- .33537943559210349312*10 + 16785.676157258815796 I], CI] 

(b) 

GF5 := [ 

CI = [ - .11869586429575443972*10 I, 

-5 -5 
.37289405527131873878*10 +.66133100616139736842*10 I], 

C2 = 

[.052942866215076463492 I, - .16632491956127147989 - .32145019412334161258 I], 
CI (1 &* ONE) + C2 (1 &* DgO)] 

(c) 

GF6 := [ 
CI = 

[ - 6.4357348171046949024 I, 20.218457221868161079 + 34.784995718164094540 I], 
C2 = [.00077427320705022066450 I, 

- .0024324510191403821401 - .0041849285213475828278 I], 
CI G(nul, mi2) + C2 ql(nul) ql(nu2)] 

Figure (a) shows a scalar self-energy. So the result just consists of the scalar constant CI. 
Like in the previous example CI is expanded as a series in e. Figure (b) shows a fermionic 
self-energy, which has Dirac structure. So the form factors are the unity matrix ONE and 
DgO. The pre-factors CI and C2 are again represented as lists of coefficients of the Laurent 
series in e. Finally figure (c) shows a vector boson self-energy. The result is proportional 
to the metric tensor G(nul, nu2) and to the product of the external momenta ql(nul) 
ql(nu2). Again two constants, CI and C2, are needed. 

A. 2. 2 Two-loop diagrams 

Now some two-loop diagrams shall be evaluated numerically. The values are chosen as 
in the previous section, if no comment is made. First we want to show some two-loop 
diagrams which factorize, starting with a photon self-energy: 
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w w 




w w 



(q10,0) (qio.o) 

>■ > 



GF7 := [CI = [ - .16744033026458920315 I, 

-19 

.70833637543010397531*10 + 2.1988161494279602884 I, 

-17 

- .42232923827218965144*10 - 14.715169445826507291 I], 
C2 = [.000020136974316008130172 I, 

-23 

- .85187071577025180931*10 - .00026443751189864126887 I, 

-21 

.50790828055269934758*10 + .0017696990248291598135 I], 
CI G(nul, nu2) + C2 ql(nul) ql(nu2)] 

The typical structure of a photon self-energy, 

CI G(nul, nu2) + C2 ql(nul) ql(nu2), 

can be seen. The difference to the one-loop self-energy (c) is just the fact, that this graph 
is more divergent. So the constants CI and C2 have an additional list entry, the order e~ 2 . 

Next we demonstrate a three-point function. All momenta are on-shell, e.g. ql0:=Mz0; 
q20:=MzO/2; q21 : =sqrt (Mz0~ 2-Melec~2) ; is passed to Maple with Insert Maple Com- 
mand. 



(q10-q20,-q21) 
>• 

e 




(q10,0) (q20,q21) 
>■ *" 



GF8 := [ 

-35 -25 
[CI = [0, - .82353974706799645789*10 - .64017338869752026115*10 I, 
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-8 -9 
- .30410069590311173109*10 - .86001332568491139242*10 I], 

-34 -25 
C2 = [0, .16428685998560205727*10 - .72770813279248241968*10 I, 

-7 -8 

- .31775075936778308993*10 - .89861644838148378070*10 I], 

-27 -27 
C3 = [0, .25907616430369099002*10 + .17084563249491568053*10 I, 

-10 -10 
.50023692425121327568*10 -.38336659269131090697*10 I], 

-27 -26 
C4 = [0, .56375581291706017052*10 + .83584350815468097846*10 I, 

-9 -9 
.52269088721608778944*10 -.40057463723369682265*10 I], 

-27 -27 
C5 = [0, .11444962108462501382*10 - .21574341953884228136*10 I, 

-26 -26 

- .13104643398124069502*10 + .25219418322926334904*10 I], 

-27 -27 
C6 = [0, .19581694212974079513*10 - .32129543956635657630*10 I, 

-26 -26 

- .22477792468575032187*10 + .37675015987648266663*10 I], 

-31 -31 
C7 = [0, .13476989938404571673*10 + .21449295663578860996*10 I, 

-23 -14 
.10998568100*10 - .26575450025692144978*10 I], 

-31 -31 
C8 = [0, - .13476989941451368157*10 - .21449295669672453982*10 I, 

-22 -13 

- .11492263298*10 + .27768333121127352587*10 I], 

-26 -26 
C9 = [0, .31582291657035101225*10 - .30101578572250827460*10 I, 

-9 -10 
.10838466684439346062*10 + .56893589844240048516*10 I], 

-36 -26 
CIO = [0, - .98315817920409517794*10 - .13705392377889147738*10 I, 

-8 -9 
.11324969215000662291*10 + .59447352866059701005*10 I], 

-27 -26 
Cll = [0, .76454675644007305666*10 - .48786422812941166168*10 I, 
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9 



-9 



.25011846212167898923*10 



- .19168329634264528256*10 



I], 



-26 



-25 



C12 = [0, - .20628807340956139895*10 



- .21423923511009696554*10 



I. 



-8 



-8 



.26134544360394034972*10 - .20028731861370353818*10 I], 
CI (1 &* Dg(nul)) + C2 &*(1, Dg5, Dg(nul)) + C3 (1 &* DgO) q2(nul) 
+ C4 &*(1, Dg5, DgO) q2(nul) + C5 (1 ft* Dgl) q2(nul) 
+ C6 &*(1, Dg5, Dgl) q2(nul) + C7 (1 ft* ONE) ql(iml) 
+ C8 (1 ft* Dg5) ql(nul) + C9 (1 &* DgO) ql(nul) 
+ CIO &*(1, Dg5, DgO) ql(nul) + Cll (1 ft* Dgl) ql(nul) 
+ C12 &*(1, Dg5, Dgl) ql(nul) ] 

Of course the result has more form factors as the diagrams evaluated before. The reason 
for this is, that as external particles vector bosons as well as fermions occur. 
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B.l File structure 




.xloops 



<home directory> 




<process name> 




Init 

History.tab 
GM.ma 



xloops_user 



Grn.ma 



.cxx 



To understand which files tfoops uses for loading and saving of processes, it is necessary 
to have a short glance on the file structure, which Troops writes into every user's home 
directory. First of all sloops creates the subdirectory xloops_user. In this directory sloops 
writes all user-specific data. If the user now starts to input a new process, xloops creates 
in xloops_user a new subdirectory, which has the same name as the process. In this 
subdirectory Troops writes to basic files: the file Init with all options entered for this process 
(e.g. the model), and the file History.tab, where for every diagram a corresponding 
number and a particle list is saved. If one graph from this list is calculated, ^loops creates 
a new file to save the result. The name of this file includes the number of the diagram 
in the particle list. If the user now tries to evaluate this graph again, Troops looks for the 
corresponding file and reads the result. 

B.2 The one-loop library 

To avoid multiple evaluation of one-loop integrals \loops can read these integrals from a 
pre-built library. This library is produced with by the OneLoopLib procedures and written 
to the directory, which is indicated by the variable LibPath. If the Xwindows interface is 
used, this variable points to a common directory, where all integrals are saved. Otherwise 
the user can set this variable to the directory, where he has his own library stored. If some 
integrals in the library are missing, sloops creates them automatically.^ 

In practice it is necessary to store not only the general mass case of each integral. 
Several special cases are also needed. The correspondence is 




^PoPiP2tit2Hnin2hi2h'iii^ ■ 1L00P <^=^ 
0neLoop3Pt (p 0j Pi , P2 , 1i , 120 , #21 , m\ , m 2 , m 3 , t\ , t 2 , t 3 ) . 



Currently only supported on Unix platforms. 
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The additional numbers 77-1, 77-2 and the optional numbers ii, 12, 13, 14, 15, i$ specify the dif- 
ferent kinematical arrangements which have to be taken into account. In the following 
tables we list the conditions which belong to each case: 

Two-point functions: 



conditions 




method of 


effective 


1 


2 


m 


solution 


n-point 


q = 


nil = TU2 


triv. 


D 


1 








P 


1 


nil = 




1 


R 


2 


ni2 = 




2 


R 


2 


q 2 + m{ — m'2 = 




3 


R 


2 


q 2 — mf + m 2 = 




4 


R 


2 


q + mi + rn 2 = 




5 


R 


2 


q — nil + = 




6 


R 


2 


q + nil — m 2 = 




7 


R 


2 


q — nil — n-2 = 




8 


R 


2 



Three-point functions: 





conditions 




method of 


effective 


1 


2 


3 


4 


77i77 2 


solution 


77-point 


921 = 


920 = 


91 =0 


"1-1 = 771-2 = 177-3 


triv. 


2 x D 


1 








"ii = m-2 


triv. 


P D 


1 








1?7i = 7173 






1 








177 2 = m-3 






1 










triv. 


P 


1 






m 2 = m 3 


gi + 7771 + 7773 = 


11 


D R 


2 








qi + m-i - rn 3 = 


12 












(/l — 777-1 + 1Tl3 = 


13 












qi — mi — 1113 = 


14 














10 


D 


2 




91 = 920 


nil = m-2 


gi + 7771 + 7773 = 


21 


D R 


2 








qi + m-i - rn 3 = 


22 












gi — mi + m3 = 


23 












qi — mi — wi3 = 


24 














20 


D 


2 




920 = 






40 


P 


2 




91 =0 






50 


P L 


2 




9i = 920 






60 


P 


2 




(920 - m 2 


91 - {qf ~ 


1^1)920 - ^3(920 - 9i) = 


70 


P 


2 










80 


P 


2 


9i =0 


mi = m-3 


Vlio ~ Q21 + m i + m 2 = 


31 


L D R 


2 






\Aio - 921 + m i -m 2 = 


32 










\Aio ~ ?2i - m i + m 2 = 


33 










\Aio ~ ?2i - m i - m 2 = 


34 














30 


L D 


2 










90 


P L 


2 
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CO 

1 


nditions 
2 


3 


4 


77i77 2 


method of 
solution 


effective 
77-point 




H2U yzl 










3 


777.3 7^ 






CO 


S 


3 


920 - 91 = -921 


920 = 921 










3 


7773 / 






DO 


S 


3 


920 = 921 


777 2 / 






EO 


S 


3 


920 = -921 


7772 






FO 


S 


3 


920 - 92i = m'i 


q\ = mf 


7773 = 




GO 


S 


3 


(920 - 9i) 2 - 92i = m i 


920 - 921 = m 3 


7772 = 




HO 


S 


3 


(920 - 9l) 2 - 921 = m 2 


'1 2 

91 = "13 


777i = 




10 


S 


3 



The table shows in the first columns the conditions to the masses rrii and the external 
momenta for a special case. Then the number occuring in the name of the file, in which 
the special case is saved is given. In the following two columns the method of solving this 
integrals and the effective 77-point function to solve is given. The abbreviations for the 
method of solution have the following meaning: 

D Differentiation to the squares of masses. 

P Solving by partial fraction decomposition. 

R Special case of the 1Z function, either vanishing or equal arguments. 
L Solution by Lorentz transformation and 

reduction of parallel space. 
S Summands vanish in the sum of residua, 

e.g. some residua are 0. 
For the tensor integrals of the three point function the following additional cases occur: 



conditions 






method of 


effective 


1 




2 


3 




solution 


77-point 


777i = 








ii = 4 


R 


3 






7772 = 




i\ = 6 


R 


3 








7773 = 


h = 7 


R 


3 






7773 = 




i\ = 5 


R 


3 


7772 = 








h = 2 


R 


3 






7773 = 




i\ = 3 


R 


3 


7773 = 






h = 1 


R 


3 


q\ - m\ + 7773 = 




777i 7^ 




i 2 = 1 


R 


3 


q\ + — 7773 = 








i 2 = 2 


R 


3 


q\ — m'2 + 7773 = 




777 2 / 




»3 = 1 


R 


3 


gf + 7772 - 7773 = 








i 3 = 2 


R 


3 


920 - m 2 + m t = 




777 2 / 




74 = 1 


R 


3 


Q20 + ml-m 2 3 = 








74 = 2 


R 


3 


920 - 921 - m 2 + m i = 









74 = 3 


R 


3 


920 - 921 + m 2 - m l = 









74 = 4 


R 


3 


q\ — m'f + 777,2 = 




777i 7^ 




75 = 1 


R 


3 


q\ + m\ — 777,2 = 








15 = 2 


R 


3 


(920 - 91 ) 2 - ml + ml 


= 






75 = 3 


R 


3 


(920 - 9i) 2 + mf - m| 


= 






75 = 4 


R 


3 


920 - 921 - m \ + m 2 = 





7772 / 




7 6 = 1 


R 


3 


920 - 92i +m\-ml = 









7 6 = 2 


R 


3 
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To make the notation clearer, we give the following examples: 



BpoPihh ■ 1L00P 
Bp piht 2 l.lL00P 
Bp piht 2 2. 1L00P 
BpoPlht 2 3. 1L00P 
BpoPl*l*24. 1L00P 

CpoP\p 2 tit 2 t 3 ■ 1LDDP 

CpoPlP2ht 2 t 3 2 ■ ildop 
CpoPlP2tit 2 h3 . 1L00P 
CpoPiK^i^^-lLOOP 
ZpoPiP2ht 2 t 3 b ■ 1LDDP 
CpoPiP2tit 2 t 3 6 ■ 1L00P 
CpoPiP2ii*2*37.1LDDP 
CpoPiP2iii2*38.1LOOP 
CpoPlP2ht 2 t 3 9 ■ 1LDDP 



0neLoop2Pt(p ,Pi, 9, mi, m 2 , ii, *2) 
0neLoop2Pt(p ,Pi, 9, 0, m 2 ,ti,t 2 ) 
□neLoop2Pt(p , pi, q, mi, 0, ti , t 2 ) 
DneLoop2Pt(pO)Pi) y m| — m|, mi,m 2 , ti,t 2 ) 
OneLoop2Pt(j>o,Pi, y mf — m|, mi,m 2 , t\,t 2 ) 



DneLoop3Pt(po,Pi,P2, 
DneLoop3Pt(p ,Pi,P2, 
DneLoop3Pt(p ,Pi,P2, 
0neLoop3Pt (p , Pi , P2 , 
DneLoop3Pt(p ,Pi,P2, 
0neLoop3Pt(p ,Pi,P2, 
DneLoop3Pt(p ,Pi,P2, 
DneLoop3Pt(p ,Pi,P2, 
DneLoop3Pt(p ,Pi,P2, 
DneLoop3Pt(p ,Pi,P2, 



qi,q20,q2i,mi,m 2 ,m 3 ,ti,t 2 ,t 3 ) 
qi,q 20 , 0, m 1 ,m 2 ,m 3 ,ti,t 2 ,t 3 ) 

91,920,921,91, y/q'io ~ 92i>°>*i,*2,*3) 
9i, 920, 9i - q20,m 1 ,m 2 ,m 3 , h,t 2 ,t 3 ) 
91,920, -q20,m 1 ,m 2 ,m 3 ,t 1 ,t 2 ,t 3 ) 
91,920,920 - qi,m 1 ,m 2 ,m 3 ,t 1 ,t 2 ,t 3 ) 
9i, 920, 920, mi,m 2 ,m 3 ,ti,t 2 ,t 3 ) 
qi,0, 0, m 1 ,m 2 ,m 3 ,t 1 ,t 2 , t 3 ) 
q 1 ,q 1 ,0,m 1 ,m 2 ,m 3 ,t 1 ,t 2 ,t 3 ) 
0,q 2 o,q 2 i,m 1 ,m 2 ,m 3 ,t 1 ,t 2 ,t 3 ) 



B.3 Run time and numerical stability 

Typical run times for the different integrals are displayed in the following table. We used 
two different systems: (i) a DEC Alpha 8400 workstation (ii) a PC with i586 chip, 32 MB 
RAM and 133 MHz frequency running Linux. All systems are working with Maple V. 
Without using the library we found the following run times: 





DEC Alpha 8400 


i586 133MHz 


tensor rank 


numerical 


algebraical 


numerical 


algebraical 


One-loop two-point functions 





0.02 s 


0.04 s 


0.04 s 


0.06 s 


1 


0.04 s 


0.21 s 


0.06 s 


0.26 s 


2 


0.07 s 


0.28 s 


0.07 s 


0.48 s 


3 


0.08 s 


0.41 s 


0.09 s 


0.76 s 


One-loop three-point functions 





2.54 s 


1.83 s 


4.83 s 


3.31 s 


1 


2.57 s 


2.23 s 


5.09 s 


3.69 s 


2 


3.01 s 


3.62 s 


5.67 s 


6.14 s 


3 


3.78 s 


4.16 s 


7.61 s 


34.58s 



It should be emphasized that these are the times which are necessary to generate the 
functions once and forever using the DneLoopLibnPt routines. If they once are stored in 
the LibPath directory they may be read in quickly. With assigned LibPath we get: 
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DEC Alpha 8400 


i586 133MHz 


tensor rank 


numerical 


algebraical 


numerical 


algebraical 


One-loop two-point functions 





0.34 s 


0.03 s 


0.08 s 


0.02 s 


1 


0.39 s 


0.05 s 


0.09 s 


0.03 s 


2 


0.43 s 


0.06 s 


0.19 s 


0.03 s 


3 


0.47 s 


0.06 s 


0.20 s 


0.04 s 


One-loop three-point functions 





1.30 s 


0.09 s 


4.77 s 


0.47 s 


1 


1.32 s 


0.10 s 


4.81 s 


0.57 s 


2 


1.36 s 


0.12 s 


5.20 s 


0.74 s 


3 


1.88 s 


0.45 s 


6.72 s 


1.71 s 



The accuracy and stability of numerical results may suffer from cancellations of large, 
approximately equal dilogarithms. Since Maple supports calculations of arbitrary precision, 
it is always possible to increase the number of digits to improve accuracy - instead of 
Fortran programs which are limited to Fortran's restricted accuracy. 

Usually we calculated with 20 or 40 digits. For all practical purposes our numerical 
results were not influenced significantly - only in the last four or five digits - if the number 
of digits increased. 



B.4 The \Coops distribution 

The distribution (packed in xloops . tgz) consists of the following files: 



README 
XLwidget .tcl 
configure* 
convert .tcl 
graphen.tcl 



graphenO.tcl 
graphenl .tcl 
graphen2.tcl 
graphen2a.tcl 
graphen2b.tcl 



graphenp.tcl 
help. tcl 
history .tcl 
maple. tcl 
model .tcl 



particles .txt 
xlinit .tcl 
xloops .prototyp 



MapleVRl : 
cf cn.ma 
evalproc .ma 
fmrules .ma 
fmuser .ma 
loops .ma 



make_lib.map 
mess . isolatin 
mess. ma 
mess.no_iso 
numint . ma 



numint .ma_unix 
numint . ma_vms 
oneloop.ma 
pv.ma 
r .ma 



simple .ma 

test_maple_bug.map 
twoloop.ma 
unvalue .ma 
values .ma 



MapleVR3 : 
cf cn.ma 
evalproc .ma 
fmrules .ma 
fmuser .ma 
loops .ma 



make_lib.map 
mess . isolatin 
mess .ma 
mess.no_iso 
numint . ma 



numint .ma_unix 
numint . ma_vms 
oneloop.ma 
pv .ma 
r .ma 



simple .ma 

test_maple_bug.map 
twoloop.ma 
unvalue .ma 
values .ma 



divO . cxx 


gauss. h 




makehead . Solaris 


numint2 . vms* 


divO.h 


incl . cxx 




makethread.O 


nvegas . c 


dpcl . in 


main . cxx 




makethread. 1 


pvegas . c 


dpc2 . in 


makefile 


rest 


makethread.2 


utils . cxx 


dpc3 . in 


makefile 


vms 


numint . cxx 


utils.h 


dpc4 . in 


makehead 


AIX 


numint .prototyp 


vegas . c 


f unct . cxx 


makehead 


Linux 


numint 1 .unix* 




funct .h 


makehead 


generic 


numint 1 .vms 




gauss . cxx 


makehead 


osf 


numint2 . unix* 
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lib: 

BOOllO.lloop 



manual : 
manual . html 



xbms : 

HpllltO.xbm 

Hpll2t0.xbm 

Hpll2tl.xbm 

Hpll2t2.xbm 

Hp210t0.xbm 

Hp211t0.xbm 

Hp211tl.xbm 

Hp212t0.xbm 

Hp212tl.xbm 

Hp212tl0.xbm 

Hp212tll.xbm 

Hp212t2.xbm 

Hp212t3.xbm 

Hp212t4.xbm 

Hp212t5.xbm 

Hp212t6.xbm 

Hp212t7.xbm 

Hp310t0.xbm 

Hp311t0.xbm 

Hp311tl.xbm 

Hp312t0.xbm 

Hp312tl.xbm 

Hp312tl0.xbm 

Hp312tll.xbm 

Hp312tl2.xbm 

Hp312tl3.xbm 

Hp312t2.xbm 

Hp312t3.xbm 

Hp312t4.xbm 

Hp312t5.xbm 

Hp312t6.xbm 

Hp312t7.xbm 

Hp312t8.xbm 



Hp312t9.xbm 

Hp410t0.xbm 

Hp411t0.xbm 

Hp411tl.xbm 

Hp411t2.xbm 

Hp412t0.xbm 

Hp412tl.xbm 

Hp412tl0.xbm 

Hp412tll.xbm 

Hp412tl2.xbm 

Hp412tl3.xbm 

Hp412tl4.xbm 

Hp412tl5.xbm 

Hp412tl6.xbm 

Hp412tl7.xbm 

Hp412tl8.xbm 

Hp412tl9.xbm 

Hp412t2.xbm 

Hp412t20.xbm 

Hp412t21.xbm 

Hp412t22.xbm 

Hp412t23.xbm 

Hp412t24.xbm 

Hp412t25.xbm 

Hp412t26.xbm 

Hp412t27.xbm 

Hp412t28.xbm 

Hp412t3.xbm 

Hp412t4.xbm 

Hp412t5.xbm 

Hp412t6.xbm 

Hp412t7.xbm 

Hp412t8.xbm 



Hp412t9 . xbm 

attention. xbm 

file .xbm 

help. xbm 

insmath . xbm 

leftarrow.xbm 

pllltO.xbm 

pllltl .xbm 

pll2t0.xbm 

pll2tl.xbm 

pll2t2.xbm 

pll2t3.xbm 

plproc .xbm 

p210t0.xbm 

p210tl.xbm 

p211t0.xbm 

p211tl.xbm 

p211t2.xbm 

p211t3.xbm 

p212t0.xbm 

p212tl.xbm 

p212tl0.xbm 

p212tll.xbm 

p212t2.xbm 

p212t3.xbm 

p212t4.xbm 

p212t5.xbm 

p212t6.xbm 

p212t7.xbm 

p212t8.xbm 

p212t9.xbm 

p2proc .xbm 

p310t0.xbm 



p310t2.xbm 

p311t0.xbm 

p311tl.xbm 

p311t4.xbm 

p311t5.xbm 

p312t0.xbm 

p312tl.xbm 

p312tl0.xbm 

p312tll .xbm 

p312tl2.xbm 

p312tl3.xbm 

p312t2.xbm 

p312t3.xbm 

p312t4.xbm 

p312t5.xbm 

p312t6.xbm 

p312t7.xbm 

p312t8.xbm 

p312t9.xbm 

p3proc .xbm 

p410t0.xbm 

p410t3.xbm 

p411t0.xbm 

p411tl.xbm 

p411t2.xbm 

p411t6.xbm 

p411t7.xbm 

p411t8.xbm 

p412t0.xbm 

p412tl.xbm 

p412tl0.xbm 

p412tll .xbm 

p412tl2.xbm 



p412tl3.xbm 

p412tl4.xbm 

p412tl5.xbm 

p412tl6.xbm 

p412tl7.xbm 

p412tl8.xbm 

p412tl9.xbm 

p412t2.xbm 

p412t20.xbm 

p412t21 .xbm 

p412t22.xbm 

p412t23.xbm 

p412t24.xbm 

p412t25.xbm 

p412t26.xbm 

p412t27.xbm 

p412t28.xbm 

p412t3.xbm 

p412t4.xbm 

p412t5.xbm 

p412t6.xbm 

p412t7.xbm 

p412t8.xbm 

p412t9.xbm 

p4proc .xbm 

rightarrow . xbm 

schraf f . xbm 

xloops .gif 

xloops .xbm 
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